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Abstract 
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iments at the front-end of a muon storage ring. We estimate the forseeable ac- 
curacy in the determination of several interesting observables, and explore the 
consequences of these measurements. We discuss the extraction of individual 
quark and antiquark densities from polarized and unpolarized deep-inelastic 
scattering. In particular we study the implications for the undertanding of the 
nucleon spin structure. We assess the determination of Os from scaling vi- 
olation of structure functions, and from sum rules, and the determination of 
sin^ 9w from elastic ue and deep-inelastic up scattering. We then consider the 
production of charmed hadrons, and the measurement of their absolute branch- 
ing ratios. We study the polarization of A baryons produced in the current and 
target fragmentation regions. Finally, we discuss the sensitivity to physics be- 
yond the Standard Model. 



CERN-TH/2001-131 
February 1, 2008 



Contents 



1 


INTRODUCTION 








2 


GENERALITIES 





3 UNPOLARIZED STRUCTURE FUNCTIONS 

3.1 Formalism 

3.2 Leading-order results 

3.3 NLO extraction of parton densities 



4 POLARIZED STRUCTURE FUNCTIONS 

4. 1 Formalism 

4.2 Positivity bounds on polarized densities 

4.3 Theoretical scenarios and the spin of the proton . . . . 

4.4 Statistical errors on polarized densities at the ;^-Factory 

4.5 Structure-function fits 

4.5.1 Results for first moments 

4.5.2 Results for x distributions! 

4.6 Extraction of As{x) and As{x) from tagged charm . . 



5 MEASUREMENTS OF as 

5. 1 Determination of ag and higher-twist terms from QCD fits of the SF date 

5.2 Determination of gg from the Gross-Llewellyn Smith sum rule 

5.3 Measurement of Fi{x) and unpolarized Bjorken sum rule 



5 NUCLEAR EFFECTS IN DIS AT THE i^-Factory 

6.1 Nuclear shadowing 

6.2 Nuclear effects at large x 

6.3 Nuclear effects in DIS sum rules 



7 ELECTROWEAK STUDIES AT THE ;y-FACTORY 

7.1 ue^ scatterin"3 

7.2 sin^ from DIS 

7.3 New physics through radiative corrections .... 



8 STUDIES WITH HEAVY QUARKS 

8.1 Direct evaluation of the A+ branching ratios 

8.1.1 Model-independent extraction of BR{A^ pK^ir^) 

8.1.2 Present kno wledge, theoretical and experimental, of neutrino quasi-elastic charm 
production 

8.1.3 Description of the method 

8.1.4 Measurement accuracy 

8.2 Direct evaluation of branching ratios and f/)^ measurement 



1 



8.2.1 Topology of neutrino-induced diffractive charm events and background 

8.2.2 Description of the method 

8.2.3 Measurement accuracy at a neutrino factor)] 

8.3 Theoretical estimates for i/-induced exclusive production 



9 A POLARIZATION IN NEUTRINO PIS 

9.1 A polarization in the current fragmentation region 

9.1.1 Charged current neutrino processes, up^iA ' X 

9.1.2 Neutral current neutrino processes, up^uA^ X 

9. 1.3 Present knowledge on A fragmentation functions and numerical results 

9.2 Models for A polarization in the target fragmentation region 

9.2.1 SC/(6) Quark-diquark model 

9.2.2 Meson-cloud model 

9.2.3 Polarized-intrinsic-strangeness model 

9.3 Discussion 



10 SEARCHES FOR NEW PHYSICS 

10.1 A search for Z' in muon-neutrino-associated charm production 

10.1.1 The process 

10.1.2 Description of the method 

10.1.3 Measurement accuracy 

10.2 Bounds on 4-fermion operators from a i/-Factory 



U CONCLUSIONS 



PHYSICS AT THE FRONT-END OF A NEUTRINO FACTORY: A 
QUANTITATIVE APPRAISAL 



M.L. Mangano " (convener), S.I. Alekhin ^, M. Anselmino R.D. Ball M. Boglione ^ U. D'Alesio f, 
S. Davidson S', G. De Lellis ^, J. Ellis S. Forte P. Gambino T. Gehrmann °, A.L. Kataev 
A. Kotzinian S.A. Kulagin ^, B. Lehmann-Dronke K P. Migliozzi ^, F. Murgia ^ , G. Ridolfi ""'"^ 
" Theoretical Physics Division, CERN, Geneva, Switzerland 

* Institute for High Energy Physics, Protvino, Russia 

Dipartimento di Fisica Teorica dell'Universita e Sezione INFN di Torino, Turin, Italy 

Dept. of Physics and Astronomy, University of Edinburgh, Scotland 

Dept. of Physics and IPPP, University of Durham, U.K. 
f Dipartimento di Fisica deH'Universita e Sezione INFN di Cagliari, Cagliari, Italy 
^ Theoretical Physics, Oxford University, U.K. 

INFN, Sezione di Napoli, Naples, Italy 

* INFN, Sezione di Roma III, Rome, Italy 

^ Institute for Nuclear Research, Academy of Sciences, Moscow, Russia 
^ JINR, Dubna, Russia 

' Institut fiir Theoretische Physik, Univeristat Regensburg, Germany 
™ INFN, Sezione di Genova, Genoa, Italy 

Abstract 

We present a quantitative appraisal of the physics potential for neutrino exper- 
iments at the front-end of a muon storage ring. We estimate the forseeable ac- 
curacy in the determination of several interesting observables, and explore the 
consequences of these measurements. We discuss the extraction of individual 
quark and antiquark densities from polarized and unpolarized deep-inelastic 
scattering. In particular we study the implications for the undertanding of the 
nucleon spin structure. We assess the determination of as from scaling vi- 
olation of structure functions, and from sum rules, and the determination of 
sin^ 9\Y from elastic ue and deep-inelastic up scattering. We then consider the 
production of charmed hadrons, and the measurement of their absolute branch- 
ing ratios. We study the polarization of A baryons produced in the current and 
target fragmentation regions. Finally, we discuss the sensitivity to physics be- 
yond the Standard Model. 

1 INTRODUCTION 

The use of intense neutrino beams as a way of exploring the deep structure of hadrons has long been 
recog nized as a big added value of a muon-collider [Q] and neutrino-factory (z^-Factory) complex. 
Recent documents have outlined with great care the areas where deep-inelastic-scattering (DIS) 

experiments operating closely downstream of the muon ring could provide significant contributions to 
our understanding of the nucleon structure. 

These studies pointed out the potential for measurements of unparalleled precision of both un- 
polarized and polarized neutrino structure functions (SFs), leading to an accurate decomposition of the 
partonic content of the nucleon in terms of individual (possibly spin-dependent) flavour densities. In 
addition to the measurements of SFs, the large rate of charm production, allowed even with muon beam 
energies as low as 50 GeV, gives an opportunity for accurate studies of the spectrum and decay prop- 
erties of charmed systems (mesonic and baryonic), as well as for an improved determination of the 
CKM matrix element Vd- Operation at muon beam energies in excess of 500 GeV would allow similar 
studies using 6-flavoured hadrons. The large neutrino fluxes will also make large-statistics u^e and UgS 



scattering experiments possible. These measurements may provide very accurate determinations of the 
weak interaction parameter sin^ 9w, complementing in terms of accuracy and systematics the current 
determinations from higher-energy measurements in decays and from DIS. 

The goal of the work performed within our Working Group was to address some of the topics pro- 
posed in ^ in a quantitative way, and carry out a concrete appraisal of the impact that measurements 
done at the z^-Factory could have on relevant observables. A first study in this direction, limited to the 
case of SFs, has recently appeared in 

In Section 2 of this document we review our notation and describe the benchmark beam and 
detector parameters used in this study. In Section 3 we discuss the determination of unpolarized SFs, and 
of their flavour decomposition, using a next-to-leading-order (NLO) global fit analysis. The importance 
of the NLO analysis is not related to the cross-section changes induced by NLO corrections, which 
are marginal when evaluating at this stage the expected event rates, but to the mixing between quark 
and gluon contributions which arise at NLO. This mixing leads to a potential loss of accuracy in the 
extraction of individual flavours. A similar NLO study is documented in Section 4 for the polarized 
case; there we study both the accuracy in the determination of the individual shapes of polarized parton 
distributions, and the accuracy in the extraction of the proton axial charges. We shall put these results in 
the framework of the ability to distinguish between different scenarios for the description of the proton 
spin. The relevance of the NLO effects is even more significant in this case than in the unpolarized case, 
because of the larger uncertainties on the polarized gluon contribution. In that section we also analyse 
the use of tagged charm final states to study the strange quark polarized distribution. In Section 5 we 
discuss the prospects for extractions of as from global SF fits, as well as from the GLS and unpolarized 
Bjorken sum rules, and in Section 6 we analyse the nuclear effects involved in the extraction of charged- 
current (CC) neutrino SFs from heavy targets. New prospects in this area are opened by the availability 
of new SFs, whose nuclear corrections have sizes different from those studied with data available today. 
In Section 7 we discuss the extraction of sin^ 9\y from ue scattering and DIS. The large statistics will 
enable measurements of an accuracy similar to that available today from LEP, and will provide important 
and complementary tests of the Standard Model (SM). In Section 8 we study measurements involving 
charm quarks. In Section 9 we consider the application of polarization measurements in semi-exclusive 
final states to the study of polarized nucleon densities. In Section 10 we finally consider the potential of 
the i/-Factory for the detection of indirect evidence of new physics from precise measurements of SM 
observables. 

2 GENERALITIES 

For our studies (and unless otherwise indicated) we shall assume the following default specifications. 
Muon beam energy, E'^ = 50 GeV; length of the straight section, L = 100 m; distance of the detector 
from the end of the straight section, d = 30 m; number of muon decays per year along the straight section, 
= 10^^; muon beam angular divergence, 0.1 x mf^/E^, being the muon mass; muon beam 
transverse size ax = (Jy = 1.2 mm. We also assume a cylindrical detector, with azimuthal symmetry 
around the beam axis, with a target of radius ii = 50 cm and a density of 100 g/cvo? (10 g/cm^ in 
the case of polarized targets). The statistics, then, scale linearly with the detector length, while the 
dependence of other parameters, such as the radius or the length of the straight section, is clearly more 
complex. Some examples are given in Fig. |l|. 

The neutrino spectra are calculated using standard expressions for the muon decays (see e.g. |^). 
For simplicity (with the exception of the ve scattering studies), we shall confine ourselves to the case 
of and CC DIS. The laboratory-frame neutrino spectra, convoluted with the CC interaction cross- 
sections, are shown for several detector and beam configurations in Fig. ^ {E^ = 50 GeV) and Fig. ^ 
{Efj^ = 100 GeV). The number of events, in different bins of (x, Q^), are shown in Fig. ^ 
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Fig. 1: Bottom: for 50 GeV muons decaying along a straight section of length L — 100, 200 and 1000m, we plot the fraction 
of the muon-neutrino flux contained within a circle of radius -R, at a distance of 30m from the end of the straight section. Top: 
average V/j, energy for the three beam configurations, as a function of R. 




K (GeV) E„ (GeV) 

Fig. 2: CC event rates, in units of 10^, as function of Lab-frame neutrino spectra, for several detector and beam configurations. 
The dashed lines on the left include cuts on the final-state muon (Ef^ > 3 GeV) and on the final-state hadronic energy (Ehad > 
1 GeV). The solid lines have no energy-threshold cuts applied. The three set of curves correpsond to different detector radiuses 
(50, 10 and 5 cm, from top to bottom). 
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Fig. 4: Event rates, in different bins of (x, Q^), for the default beam and detector configuration (Ey, — 50 GeV, L = 100 m, 
i? = 50 cm). 
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Fig. 5: X and binning for the generation of CC events. The crosses correspond to the weighted bin centers. 



3 UNPOLARIZED STRUCTURE FUNCTIONS 
3.1 Formalism 

Unpolarized CC SFs are defined through the decomposition of unpolarized differential CC cross-sections 
into invariant functions of the momentum of the struck quark (x) and the momentum transfer squared of 
the W boson {Q^): the standard definitions give 



dxdy 27r(l + QVMj 



W) 



(1 - y)F2 + y'xFi ± y ( 1 - | ) xFg 



(1) 



where S = 2mEy is the nucleon-neutrino centre-of-mass energy, m is the nucleon mass, is the 
neutrino beam energy, assumed to be S> m, y is the fractional lepton energy loss, or {By — E£)/Ey, 
and the it signs refer to the sign of the CC: exchange for v scattering and for u. In neutrino 
scattering, x, y, and can all be determined simply by measuring the outgoing lepton energy and 
direction, and the hadronic energy in the event. If both v and v beams are available, there are then, in 
principle, six independent SFs for each target. 

We now wish to determine the expected statistical accuracy with which the individual SFs, and 
their flavour components, can be determined. To do this, we shall exploit the different y dependences 
of the cross section on the various Fj. The advantage of the neutrino beams from muon decays is their 
wide-band nature. This allows us to modulate the y dependence for fixed values of x and using the 
neutrino energy: 

... 

We produced y distributions by generating events within different bins of x and Q^, and performed 
minimum-x^ fits of the generated data using the cross-section Eq. (|I]). For each bin, the values of x and 
at which we quote the results are obtained from the weighted average of the event rate. As an input, 
we used the CTEQ4D set of parton distributions [^. The dependence on the parameterization of the 
parton distributions is very small, and will be neglected here. We verified that other recent sets of parton 
distributions give similar results. The absolute number of events expected in each bin is scaled by the 
total number of muon decays; this number of events determines the statistical error on the individual SFs 
obtained through the fit. 

We generate events in the (x, Q"^) bins shown in Fig. |5[ Twenty equally-spaced bins in the range 
< y < 1 are used for the y fit. The total number of x bins varies in different bins because of 
kinematic acceptance and minimum energy cuts. The statistical errors returned by the fits are used as 
estimates of the statistical errors in the extraction of the Fj. In the parton model, four of the SFs are 
related through the Callan-Gross relations F2 = 2xFi: the longitudinal SF F^ = F2 — 2xFi begins at 
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Fig. 6: Expected errors on the determination of the individual nucleon SFs {x, Q ) (i = 1, 2, 3), for various ranges of Q . 
The horizontal bars indicate the range in x defining the bins within which the statistical errors are determined. 

0{as) in perturbation theory. We considered both the cases of three-component fits (leaving Fi and F2 
uncorrelated, and fitting all three SFs), and two-component fits (assuming the Callan-Gross relation, and 
fitting for F2 and F3). 

The three-component fit for the nucleon SFs, obtained assuming one year exposure of a deuterium 
target to a muon-neutrino beam, is shown in Fig. |[ Notice that the errors of Fi are always better than 
10%. F2 is determined very well at large x since in this region one is more sensitive to the y—>^0 limit, 
where the contribution from Fi and F3 is suppressed (see Eq. [I]). Figure ^ shows the result obtained 
assuming the Callan-Gross relation. Now the relative errors are better than 1% for most regions of x and 
Q^. The significant improvement in the determination of F3 at large x results from the Callan-Gross 
constraint on Fi, and the fact that, as pointed out above, F2 is very well determined at large x. 



3.2 Leading-order results 

The parton content of the SFs Fi, F2 and F3 depends on the charge of the exchanged gauge boson. At 
leading order, we have 

F^^ = u + d + s + c, F^~ =u + d + s + c, 

F^^ = 2x{u + d+s + c), F^~ =2x{u + d + s + c), (3) 

xF^^ = 2x{-u + d + s-c), xF^~ = 2x{u -d-s + c). 

The corresponding expressions for a neutron target are obtained (assuming isospin invariance) by replac- 
ing u d. It is thus not difficult to see that, constructing appropriate linear combinations of the eight 
independent SFs (-Fj^ )p,n and (^^3'^ )p^n, it is possible to disentangle u±u, d±d and s±s, provided 
only that c± c can be determined independently, either theoretically or empirically. More explicitly, we 
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have 



-w-\ 

3 >P 

^w++w-\ 



W++W- 



+ {-Fr- 



w- 



{xF. 



where F- 



(Fr-^-),-(F, 
{Fr-''')p+{F^ 

W+±W- — liT?W+ j^pW 



W++W- 



W+-W- 



-w- 



x{u + u + d + d + s + s + c + c), 

-2x{u + u- {d + d)), 

2x{s + s-{c + c)), 

x{u — u + d — d+s — s + c — c), 

—2x{u — u — {d — d)), 

2x{s — s — (c — c)), 



F. 



W++W- 



. The first and fourth of these equations are the SFs 



W++W' 



(4) 
(5) 
(6) 
(V) 
(8) 
(9) 

and 



normally measured in neutrino scattering (though on heavy targets). The second and third 



equations allow flavour decomposition of the total q + q distributions, and the fifth and sixth equations 
allow a similar decomposition for the valence distributions. 

A detailed study of the strange sea is especially important, since this distribution is very poorly 
known at present. In a leading-order analysis one can extract the individual s{x) and s{x) at the //-Factory 
using Eqs. (Q) and (^). In order to disentangle the strange and charm contributions, it is necessary either 
to tag charm in the final state, or to assume that the charm contribution is generated dynamically by 
perturbative evolution. Assuming that the charm contribution can be either neglected or independently 
determined, a 2-year running on a nucleon target (1 year each for fi^ and /i~ beams) would result in 
errors in the determination of s{x) + s{x) as shown in Fig. ^. There we compare the size of the predicted 
uncertainties with the values of the s{x) + s{x) densities currently estimated in the analysis of the CTEQ 
group (set CTEQ4, and in the recent study by Barone, Pascaud and Zomer (BPZ, [Q]). This last 
study, based on a fit to CDHS neutrino data, finds evidence for an intrinsic strange component of the 
proton, which results in an enhancement at large x relative to the CTEQ fits. As shown in the figure, this 
evidence could be firmly established, and its size very accurately determined, using neutrino factory data. 
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Equation can then be used to estimate the errors on the determination of the difference s{x) — s{x). 
These are shown in Fig. ^, superimposed to the recent fits by BPZ. 



3.3 NLO extraction of parton densities 

The leading-order study presented in the previous subsection already provides a useful set of benchmark 
accuracies that can be achieved at the z/-Factory. However, an estimate of the precision in the extraction of 
the individual parton distribution functions (PDFs) requires a full NLO analysis. In order to estimate the 
statistical uncertainties in the extraction of PDFs from i^-Factory data, we used the errors calculated in the 
previous section to generate 8 sets of 'fake' data for the SFs ^2^3 for neutrino/antineutrino beams and for 
hydrogen/deuterium targets. The central values of the fake data were obtained using the PDFs extracted 
from the existing charged leptons DIS data [||]. The SFs i<2,3 calculated from the PDFs, parametrized as 

xdv{x,Qo) = A'^x'^'^il-xf'', xds{x,Qo) = Ajx^^-'il-xf^^ , (10) 
xuvix,Qo) = A]^x^-{l-x)'-il + j^x), ra5(x,go) = ^fx"-(l-x)''- , (11) 
xG{x,Qo) = AGx'''3{l-x)^'3{l + -ffV^ + j^x), xs{x,Qo) = Asx'''^{l-x)'''^ , (12) 

at Qq = 9 GeV^ and evolved within the NLO QCD approximation, were then fitted to these fake 
data, varying the PDFs parameters. The form of Eq. (p^ was motivated in Ref. [^, but, contrary to 
that analysis, the PDFs parameters asu, Ag, ass, bss were also released in the fit, while the parameters 
, A}^, Ag were constrained by conservation of momentum and fermionic numbers, as in Ref.[||]. The 
indices V and S are used to denote valence and sea distributions, respectively; the functions n, d, s, G 
give U-, d-, s-quarks, and gluon distributions; and s = s is assumed at this stage of the analysis. 

The statistical errors on PDFs obtained in the fit to fake data are given in Fig. [l^. For comparison, 
we present in the same figure the errors on the PDFs obtained from the analysis of Ref. [^. One can 
see that the precision in the determination of the gluon distribution has improved by about an order of 
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Fig. 10: The relative statistical errors on PDFs (%) accessible at the //-Factory (full lines). The dashed lines give statistical 
errors on the PDF set A99 obtained from the fit to the existing charged leptons DIS data ||8||. 

magnitude, and even more for the valence u- and d-quark distributions. This improvement may be ex- 
tremely helpful, for example in expanding the capabilities of the LHC in searches for new phenomena 
characterized by particles of large mass. The precision of the determination of the sea quark distribu- 
tions attainable at the z/-Factory cannot be directly compared with the corresponding errors on the PDFs 
given in Ref . ^ ; there, several parameters describing the sea quark distributions were fixed at reasonable 
values, since available data do not allow their independent extraction. In particular, the strange sea was 
assumed equal to about a half of the non-strange sea, while the behavior of the non-strange sea-quark 
distribution at small x was assumed to be universal. Thanks to the availability of eight independent com- 
binations of parton distributions, a complete separation of individual flavor distributions is now possible 
without additional constraints. This is crucial for precise studies of the flavour content of the nucleon. 
The correlation matrix for the PDFs extracted from the QCD fit to the fake z/-Factory data is given in 



Fig. yjj. Indeed, one can see that in general the absolute values of the correlation coefficients do not 
exceed 0.3 in the whole range of x. 

We now turn to a study of the potential of the i/-Factory data to determine the asymmetry of the 
strange sea, which is not accessible in neutral-current (NC) DIS experiments. In this analysis, the strange 
sea was chosen to be of the form 



Xs{x,Qo) = AsX"''' {I - X)''''' {I + CssX'^"'), Xsix,Qo) = AsX"'' {1 - xf""' {1 + CsasX" 



(13) 



which is motivated by the results of Ref. We now assume that the s — s difference is as given in 
Ref. [^, by choosing similar values of the parameters, namely 



As = 0.06, bss = 5.6, bs 



5.4, 



11000, dss = 12, ds 



7.4. 



(14) 



The parameter Uss was set equal to agd, and the parameter Csas was calculated from the constraint 



dx [s{x) — s{x)] = 0. 



(15) 
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Fig. 11: The x-dependence of the correlation coefficients p for the PDFs extracted from the analysis of the generated v- 
Factory data. The labels in the diagonal mark the rows and columns of the matrix. 

We cannot simply combine these s and s distributions with those of Ref. [S], because the xs distribu- 



tion of Eq. ( |13| ) at large x is comparable to the contribution from valence quarks, and correlated to it. 



Therefore, we repeated the fit of Ref. [^] with the s and s distributions of Eqs. ( ^3| , |14[ ). We then used 
the results of this fit to generate a new set of fake data, and repeated a QCD fit to these fake data using 
the PDFs of Eq. (12), but with the s and s distributions of Eq. (p3|). The errors on x(s — s) obtained in 



this way are given in Fig. |12|. For comparison, we recall that the error on x{s — s) obtained in the BPZ 
analysis are 0(0.002) at x ~ 0.7 (see Fig. 10); this means that an improvement of more than an order of 
magnitude in the precision may be achieved at the z/-Factory . 

Alternatively, one can choose different forms for the x{s — '§) difference; for example, it may be 
constructed using Regge phenomenology considerations. According to this approach, the behavior of 
non-singlet parton distributions at small x is governed by meson trajectories, and is generally ~ ^/x. We 
constructed two variants of the parameterization of x{s — s) based on these arguments, namely 

x{s -s){x,Qo) = AAsX^'^^il - xf'^^il + CAsx) (16) 

and 

x{s -s){x,Qo) = AAsX^'^'il - xf^'{l + CAsX + dAsX^)- (IV) 

The starting PDF set of Eq. |l^ was then modified by substituting xs ^ xs + x{s — s)/2 and xs —>■ 
xs—x{s — s) /2. The initial values of the parameters used to generate fake data were chosen as uas = 0.5, 
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Aas = 0.02 and b^s = 7.5. The initial value of the parameter d-As of Eq. (|17|) was set equal to 100 and, 
in both cases, the parameter cas was defined from the constraint of Eq. ([T5|). 

The errors on x{s — s) obtained using a Regge-like form of x{s — s) are given in Fig. One 
can see that they are quite different from those obtained by assuming the BPZ form for x{s — s). The 
main difference between Regge-like and BPZ forms, is that the strange sea for the latter is very large at 
high X, even larger than the d-quark sea. As a result, the precision of the strange-sea determination for 
the BPZ set is comparable to the precision in the determination of the valence quark, namely 0(0.1%). 
Since s » s in the BPZ fit, A(s — s)/{s — s) ^ As/s. Thus, the high precision obtained using the BPZ 
form simply originates from the fact that s is very large at high x. 

The precision of the determination of x{s — s) accessible at the i/-Factory will however always 
be at the level of 1%, or better, in the region of its maximum, regardless of the functional form used to 
parameterize this difference. The accuracy relative to the total amount of strange sea is given, for the 
different strange parameterizations, in fig. 13. 
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Fig. 12: The error bands for x{s — s) difference accessible at the //-Factory for the various forms used for the parameterization 
of this difference (full line: the BPZ-like form, dashed lines: the Regge-like form with linear polynomial factor, dotted line: the 
Regge-like form with quadratic polynomial factor). The la bands are given for the the Regge-like parameterizations and 10a 
bands for the BPZ-like one. 



4 POLARIZED STRUCTURE FUNCTIONS 
4.1 Formalism 

Polarized SFs may be defined in analogy to the unpolarized ones through asymmetries in the polarized 
cross sections (see [Q] for a recent review). The polarized cross-section difference (with a proton helicity 
Ap = ±1), 

Aa = a{Xp = -l)-a{Xp = +l), (18) 
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Fig. 13: Accuracy in tiie determination of (s — s) relative to tiie total amount of strange sea, for the different strange parame- 
terizations. 



is given by 

dxdy 



7r(l + QVm^)2 xy 

2 



m 



Xix^y'^gi + Xi2x^yg2 + | 1 - y - x^y^ '-^ \ xgs (19) 



-A^ y{2 - y)xgi - (1 - y)g4 - y^xg^ 

n2\ 



„2„.2 



m 



^ i ^ ~ 2^ ~ ^^^^^ ] ~ ^^2/^5^5 



m 



2„,2, 



where is the lepton heUcity. With these definitions [^] [|, (72 and g^ drop out in the high energy limit 
E m, and we are left with an expression of the same form as the unpolarized decomposition ([l|), but 
with — 55, F2— > — gi and 2gi. Thus we again have only three partonic SFs for each v and 



d'^Aa^'{x,y,Q^ 
dxdy 



g'p 



7r(l + QVm^)2 xy 



-Xi y{2 - y)xgi - (1 - y)gi - y'^xg^ 



(20) 



The two remaining SFs ^2 and g^ have no simple partonic interpretation and are contaminated by twist-3 
contributions: their twist-2 components are fixed by the Wandzura-Wilczek relation (giving g2 in terms 
of gi) and a similar relation gives 53 in terms of g^. They are determined by measuring asymmetries 
with a transversely polarized target. In the parton model, g^ and g^ are related by an analogue [12] 
of the Callan-Gross relation: g^ = 2xg^{\ + 0{as))- Even though this relation is violated beyond 
leading order, the SFs g^ and g^ still measure the same combination of parton distributions, albeit with 
different coefficient functions. Therefore at leading twist there are only two independent polarized SFs, 
conventionally taken to be gi and g^. 

'There are many variants in the literature: see [ |ll| | for a compilation. In particular, the conventions used here are the same 
as those used in ref. except for the signs of 54 and 35. 
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The flavour decomposition of the SFs gi and may be expressed in terms of parton densities [13] 

as 

gY^ = Au + Ad + As + Ac, g^' = Au + Ad + As + Ac, (21) 

g^^ = An - Ad - As + Ac, g^' = -Au + Ad + As - Ac , (22) 

in precise analogy with the unpolarized case. Again, by constructing appropriate linear combinations of 
all eight independent SFs (conventionally taken as {g^ )p^n and {gY^)p^n), obtained by longitudinally 
polarized u and u scattering on proton and neutron (or deuteron) targets, it is possible to separately 
disentangle Au it Au, Ad it Ad and (As it As) just as in Eqs.(^-(^, provided only that Ac ± Ac can 
be determined. 

Some combinations of polarized SFs are of particular interest. For example, writing gj^'^^^ = 
gW+ _|_ gW ^ jj^g moment of 

AS = (5r^+^~)p = {gV^^^~)n = Au + Au + Ad + Ad + As + As + Ac + Ac (23) 

is the singlet axial charge uq. This is a much more direct measurement than the traditional one through 
electron-proton or deuteron DIS, since in the latter case one must first subtract the octet charge as which 
is then only determined indirectly through hyperon decays. Thus in v-DlS one would have a direct check 
on the anomalous suppression of gq. Similarly, first moments of 



(gT )p - {aT )n\ = {97*-'^~)v - {9r-'^~)n = An + An - (Ad + Ad) = Aq^ (24) 

{9r)v + (^r )n] - \ (^r^^")^ = (55^"-^' )p + {9r-'^')n = -(as + as-) + (ac + ac) , 

(25) 



give direct measurements of the axial charge 03 and of the contribution of strange quarks to the nucleon 
spin, as would the tagging of charm in the final state. The first moment of the combination 

Ags = An + An + Ad+ Ad- 2(As + As) (26) 

is the octet axial charge ag, currently determined by hyperon decays using SU(3) symmetry, thereby 
allowing a test of SU(3) symmetry violation, and specifically a distinction between different models for 
it [21. 

Flipping the signs, we can also determine the contribution of valence quarks to the spin, since 

^^w++w-^^^^gW++w--^^ = -An + An- Ad + Ad- As + As- Ac + Ac, 
(5r-'^")p-(5r-'^~)n = -2(An-An-(Ad-Ad")), (27) 
(5r-^")p + (5r-'^')n = 2(As-As--(Ac-Ac)), 

so one could even check for intrinsic strange polarization As — As. None of these valence polarizations 
can be cleanly measured in current polarization experiments. 

It should be pointed out that the flavour separations shown above receive important contributions 
from NLO corrections, most notably from contributions from the polarized gluon density. In particular, 
gi is given at NLO by 

Q2) = AQ ^ 2[nj/2] ACg ® Ag , (28) 

where the ACj are appropriate coefficient functions. Also, beyond leading order the relation between 
parton distributions and SFs is ambiguous because of the factorization scheme ambiguity. This problem 
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is particularly relevant in the case of Eq. (|28|), because, as is by now well known ||9|], the scheme depen- 
dence of the first moment of the singlet polarized quark distributions is unsuppressed as Os^O, due to 
the fact that at leading order the first moment Ag{l) of the gluon distribution evolves as l/os- A conse- 
quence of this is that it is possible to choose the factorization scheme in such a way that the first moment 
of the singlet quark distribution is scale independent (to all perturbative orders), even though this is not 



the case in the MS scheme. This is the case e.g. in the so-called Adler-Bardeen (AB) scheme [|18|]. The 
first moment of any quark distribution in the MS and AB schemes are related by 

Agp(l, Q') = Aqt^'il) - ^Ag{l, Q% (29) 

where Ag(l,Q^) = Jq dxAq{x,Q'^). Theoretical motivations for this choice will be discussed in 
Sect. p3| below. 

It follows that care should be taken when comparing AB and MS scheme polarized parton dis- 
tributions, since the quark distributions will differ by an 0(1) gluon contribution. Furthermore, in a 
generic scheme, a meaningful determination of polarized parton distributions requires at least the inclu- 
sion of NLO corrections. The full NLO analysis of parton distributions in neutrino DIS, using the known 
anomalous dimensions [ |T6| ] and coefficient functions [p^], is presented in Ref. [10]. 



4,2 Positivity bounds on polarized densities 

Cuurently available NC polarized DIS data can be used, in conjunction with specific hypotheses on the 
form of Aq — Aq, to explore potential scenarios to be probed with polarized neutrino scattering at the 



z/-Factory. To start with, it is interesting to study the constraints set by positivity [19] on the values of 
Aq — Aq for individual flavours. At leading order in QCD, the positivity of both left- and right-handed 
quark densities implies the following obvious relations: 

\Aq{x)\ < q{x) , \Aq{x)\ < q{x) . (30) 

These constraints are shown in Fig. ^ for x{Aq — Aq), with q = u,d. The allowed region is confined 
between the two continuous lines. Use was made of the most recent polarized fits from ABFR [^] and of 
the CTEQ5 [^] unpolarized densities. The figures show that, while the assumption Aq = is consistent 
with the positivity bounds, the hypothesis Aq = Aq badly violates them as soon as x > 0.1-0.2. As a 
result we conclude that | Ag| S> \Aq\ for x > 0.1, which is reasonable in view of Eq. (|^ and of the fact 
that g » g in this region of x. 

In the case of the strange quark, it turns out that the combination As + As from the ABFR fit 
already violates the positivity bound obtained using the unpolarized strange distributions from CTEQ5. 
This bound is instead satisfied if the unpolarized strange distribution from the BPZ fit is used. In such 
a case, the bound turns out to be satisfied for both Aq = and Aq = Aq (see Fig. |^). The NLO 
corrections to the positivity bounds turn out to be negligible (see Ref. [p^ for a detailed discussion). 



4.3 Theoretical scenarios and the spin of the proton 

One of the main reasons of interest in polarized quark distributions is the unexpected smallness of the 
nucleon axial charge, which has been determined in the first generation of polarized DIS experiments. A 
clarification of the physics behind this requires a determination of the detailed polarized parton content 
of the nucleon. It is useful to sketch here the various scenarios to describe the polarized content of 
the nucleon, which are representative of possible theoretical alternatives and could be tested in future 
experiments. 

Firstly, it should be noticed that even though current data give a value of the axial charge which 



is compatible with zero, they cannot exclude a value as large as ao(10 GeV ) = 0.3 [p2|]. Also, the 



current value is obtained by using information from hyperon /? decays and SU(3) symmetry. Clearly, the 
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Fig. 14: Positivity bounds on a;(Ag — Aq) (q = u, d, s), compared with the Aq = and Ag = Ag hypotheses. Unpolarized 
densities taken from the CTEQ5 parametrization. 
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Fig. 15: Positivity bounds on x{As — As), compared with the As = and As = As hypotheses. s(x) unpolarized densities 
taken from the BPZ parametrization. 
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theoretical implications of an exact zero are quite different from those of a value that is just smaller than 
expected in quark models. It is thus important to have a direct determination of the axial charge. If a 
small value is confirmed, it could be understood as the consequence of a cancellation [El] between a large 



value of the scale-independent first moment of the quark (discussed at the end of Sect. [4. l\ ) and a large 
first moment of the gluon. In this ('anomaly') scenario the up, down and strange polarized distributions 
in the AB -scheme aie close to their expected quark-model values, so in particular the strange distribution 



is much smaller than the up and down distributions. In Ref. [|24|], this cancellation of quark and gluon 
components has been derived from the topological properties of the QCD vacuum (and thus further 
predicted to be a universal property of all hadrons). 

If instead the polarized gluon distribution is small, the smallness of the singlet axial charge can 
only be explained with a large and negative strange distribution. In this case, the scale-independent first 
moment of the singlet quark distribution is also small. This scale-independent suppression of the axial 
charge might be explained by invoking non-perturbative mechanisms based on an instanton-like vacuum 
configuration p5[]. In this 'instanton' scenario the strange polarized distribution is large and equal to the 
antistrange distribution, since gluon-induced contributions must come in quark-antiquark pairs. 

Another scenario is possible, where the smallness of the singlet axial charge is due to intrinsic 
strangeness, i.e. the C-even strange combination is large, but the sizes of As and As differ significantly 
from each other. Specifically, it has been suggested that while the strange distribution (and specifically its 
first moment) is large, the antistrange distribution is much smaller, and does not significantly contribute 



to the nucleon axial charge [26]. This way of understanding the nucleon spin structure is compatible with 



Skyrme models of the nucleon, and thus we will refer to this as a 'skyrmion' scenario [|27|]. 

Therefore, the main qualitative issues that are relevant to the nucleon spin structure are to assess 
how small the axial charge is, to determine whether the polarized gluon distribution is large, and then 
whether the strange polarized distribution is large, and whether the strange polarized quark and antiquark 
distributions are equal to each other or not. More detailed scenarios might then be considered, once the 
individual quark and antiquark distributions have all been accurately determined. For instance, while the 
up and down antiquark distributions are small, they need not be zero, and in fact they could be different 
from one another [^], just like their unpolarized counterparts appeal" to be. Investigating these issues 
could shed further light on the detailed structure of polarized nucleons. 

4.4 Statistical errors on polarized densities at the z/-Factory 

The fit to the y distributions at fixed x and for a fully polarized target gives the value of the combina- 
tions F2 ± 2xg^ and F3 ib 2gi . Polarization asymmetries are extracted by combining data sets obtained 
using targets with different orientations of the polarization. The statistical accuracies with which the 
combinations can be performed depend on the statistical content of each individual data set. Since the 
polarization asymmetries are small with respect to the unpolarized cross sections, the absolute statistical 
uncertainties on the extraction of polarized SFs will have a very mild dependence on the value of the 
polarized SFs themselves; they will be mostly determined by the value of the unpolarized SFs (which to 
first approximation fix the overall event rate), and by the polarization properties of the target. 

Therefore, for simplicity, we directly use the expected statistical errors (Tf2,F3 obtained in Sec- 



tion for the extraction of F2 and F3 form unpolarized targets, assuming in this case a target thickness 
of lOg/cm^. We then relate these to the errors on the polarized cross sections by using the following 
relation given in ^ : 

= Fl$V2 aij ^ , (31) 

where aij = 1 for = (1, 3) and aij = 1/x for = (5, 2), and where F*^p is a correction factor 
(always larger than 1) that accounts for the ratio of the target densities to H2 or D2, for the incomplete 
target polarization, and for the dilution factor of the target, namely the u (or u) cross-section weighted 
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ratio of the polarized nucleon to total nucleon content of the target. The factor of \/2 in the numerator 
reflects the need to subtract the measurements with opposite target polarization. 

The values of the uncertainties in the determination of the eight CC SFs {gi and with the 
two available beams and targets) are assigned at the cross-section weighted bin centres. To obtain the 



absolute errors on the SFs for proton and deuterium targets we use the p-butanol and D-butanol target [29] 
correction factors given in Ref. [g, namely PP = 2.6, = 1.6 and F^^ = 4.4 (for a more complete 
discussion of polarized targets and their complementary properties, see ^ and references therein). We 
have assumed a luminosity of 10^" muons decaying in the straight section of the muon ring for each 
charge, for each target, and for each polarization. Assuming that only one polarization and one target can 
run at the same time, this means eight years of run. While the number of muons may not be dramatically 
increased, the integration time can be reduced by a large factor if the target thickness can be increased 
over the conservative 10 g/cm^ assumed here, or if different targets can be run simultaneously. 

4.5 Structure-function fits 

We can now study how CC DIS data may be used to determine the polarized parton content of the 



nucleon. This study has been performed in Ref. [10|, which we summarize here. First, we need to 



make assumptions on the expected flavour content of the nucleon, implementing the theoretical scenarios 
summarized in Sect. 



We define two sets of C-even polarized densities consistent with existing NC data (one set cor- 
responding to a polarized gluon consistent in size with the 'anomaly' scenario, the second set corre- 
sponding to a vanishing gluon polarization at Qq = 1 GeV), and then to define three possible inputs 
for the C-odd densities, consistent with the expectations of the three scenarios. We shall then generate 
data according to these three alternatives, with the errors defined as in the previous section, and study the 
accuracy with which their parameters can be measured at the z^-Factory. 

We start by describing the parametrization of the C-even densities, for which we adopt the type-A 



fit of Ref. [p2|], defined as follows. The quark distributions AS+ (23), Aq^, ( P4| ) and Agg (26), and the 



polarized gluon distribution Ag at the initial scale Qq = 1 GeV^ are all taken to be of the form 

Afix,Ql)=MfVfx''fil-xffil + jfx'f) , (32) 

where the factor A// is such that the parameter r/j is the first moment of Af at the initial scale. The 
non-singlet quark distributions Aq^ and Aqg, are assumed to have the same x dependence, while the 
parameter r/g, corresponding to the first moment of Agg, is fixed to the value 7?g = 0.579 from octet 
baryon decay rates using SU(3) symmetry. Furthermore, 7s = 7g = 10, 63 = ds = 0.75, 6y: = = 1. 



All other parameters in Eq. (32) are determined by the fitting procedure. Recent data [30] from the 



E155 collaboration, and the final data set from the SMC collaboration have been added to those used in 



Ref. []22|], leading to a total of 176 NC data points. 

The best-fit values of the first moments of the C-even parton distributions at the initial low scale 
of Qo = 1 GeV^ thus obtained are listed in the first column of Table |l], together with the errors from 
the fitting procedure. The subsequent three rows of the table give the values of the first moments of 
Aq{x, Q"^) + Ag(x, Q"^) at = 1 GeV^ for up, down, and strange, obtained by combining the singlet 
and non-singlet quark first moments above. Finally, we give in the last row the value of the singlet axial 
charge oq at the scale Q"^ = 10 GeV^. 

Because the first moment of the gluon distribution in this fit is quite large, we can take this global 
fit as representative of the 'anomaly' scenario, even though the strange distribution is not quite zero. In 
order to construct parton distributions corresponding to the other scenarios we have also repeated this fit 
with the gluon distribution forced to vanish at the initial scale. This possibility is in fact disfavoured by 
several standard deviations; however, once theoretical uncertainties are taken into account a vanishing 



gluon distribution can only be excluded at about two standard deviations [22], and thus this possibility 
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Table 1: Best-fit values of the first moments for the data and pseudodata fits discussed in the text. 



Par. 


Generic fit 


= Ofit 


'Anomaly' refit 


'Instanton' refit 


'Skyrmion' refit 


VT, 


0.38 ± 0.03 


0.31 ±0.01 


0.39 ±0.01 


0.321 ±0.006 


0.324 ± 0.008 




0.79 ±0.19 





0.86 ±0.10 


0.20 ±0.06 


0.24 ±0.08 


m 


1.110 ±0.043 


1.039 ±0.029 


1.097 ±0.006 


1.052 ±0.013 


1.066 ±0.014 


V8 


0.579 


0.579 


0.557 ±0.011 


0.572 ±0.013 


0.580 ± 0.012 


riu 


0.777 


0.719 


0.764 ± 0.006 


0.722 ±0.010 


0.728 ± 0.009 


Vd 


-0.333 


-0.321 


-0.320 ± 0.008 


-0.320 ± 0.009 


-0.325 ± 0.009 


Vs 


-0.067 


-0.090 


-0.075 ± 0.008 


-0.007 ±0.007 


-0.106 ±0.008 


ao 


0.183 ±0.030 


0.284 ± 0.012 


0.183 ±0.013 


0.255 ±0.006 


0.250 ± 0.007 



cannot be ruled out on the basis of present data. The results of this fit for the various first moments are 
displayed in the second column of Table 

We can now use these parton distributions to construct the unknown C-odd parton distributions. 
We construct three sets of parton distributions, corresponding to the three scenarios of Section In 
all cases, we assume Au{x) = Ad{x) = 0. Furthermore, as the 'anomaly' set we take the 'generic' 
fit of Table |l| with the assumption As{x) = 0, the strange distribution for this set being relatively small 
anyway. As 'instanton' and 'skyrmion' parton sets we take the Ag = fit of Table [I[ with As = As 
in the former case, and As = in the latter case. The charm distribution is assumed to vanish below 
threshold, and to be generated dynamically by perturbative evolution above threshold. With these choices 
all quark and antiquark distributions are fixed, and thus all SFs can be computed. 

We generate for each of these three scenarios a set of pseudo data, by assuming the availability of 
neutrino and antineutrino beams, and proton and deuteron targets, in the (x, Q"^) bins of Fig. |5[ The data 
are gaussianly distributed about the values of the SFs at each data point in the three scenarios. We obtain 
in this way approximately 70 data points for each of the eight CC SFs. 

We proceed to fit a global set of data, which includes the original NC data as well as the generated 
CC data. We assign to the generated data the estimated statistical errors, and fit including statistical 
errors only. The errors assigned to the NC data are instead obtained, as in our original fits, by adding in 
quadrature the statistical and systematic errors given by the various experimental groups. 

The fits are performed by adopting the same functional form and parameters as in the original fit for 
the C-even parton distributions, except that the normalization of the octet C-even distribution r/g is now 
also fitted. For the C-odd parton distributions, we add six new parameters, namely the normalizations of 
the up, down and strange C-odd distributions, and three small-x exponents a (corresponding to an x" 
small-x behaviour). The shape is otherwise taken to be the same as that of the C-even quark distributions. 

4.5.1 Results for first moments 

The best-fit values of all the normalization parameters are shown in the last three columns of Table [T], 
where the rows labelled tj^ and rjs now give the best-fit values and errors on the first moments of 
Aq^ = Aq — Aq. A comparison of these values with those of our original fits leads to an assessment of 
the impact of CC data on our knowledge of the polarized parton content of the nucleon. 

First, we see that the improvement in the determination of the polarized gluon distribution is 
small, though significant. This is because the gluon distribution is determined by scaling violations, and 
the available range of at a E'^ = 50 GeV i/-Factory is limited. 

Let us now consider the C-even quark distributions. The error on the first moment of the singlet 
quark AS is reduced by the CC data by a factor of 3-5 relative to the available NC data. This improve- 
ment is especially significant since the determination of r/s no longer requires knowledge of the SU(3) 
octet component, unlike that from NC DIS, and it is thus not affected by the corresponding theoretical 
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uncertainty. With this accuracy, it is possible to experimentally refute or confirm the anomaly scenario by 
testing the size of the scale-independent singlet quark first moment. Correspondingly, the improvement 
in knowledge of the gluon first moment, although modest, is sufficient to distinguish between the two 
scenarios. 

The determination of the singlet axial charge is improved by an amount comparable to the im- 
provement in the determination of the singlet quark first moment. Its vanishing could thus be established 
at the level of a few per cent. The determination of the isotriplet axial charge is also significantly im- 
proved: the improvement is comparable to that on the singlet quark; it is due to the availability of the 



triplet combination of CC SFs given in Eq. ([24[). This would allow an extremely precise test of the 
Bjorken sum rule, and accordingly a very precise determination of the strong coupling. Finally, the octet 
C-even component is now also determined, with an uncertainty of a few per cent. Therefore, the strange 
C-even component can be determined with an accuracy better than 10%. Comparing this direct determi- 
nation of the octet axial charge to the value obtained from baryon decays would allow a test of different 
existing models of SU(3) violation [p^]. 

Coming now to the hitherto unknown C-odd quark distributions, we see that the up and down 
C-odd components can be determined at the level of few per cent. This accuracy is just sufficient to 
establish whether the up and down antiquark distributions, which are constrained by positivity to be 
quite small, differ from zero, and whether or not they are equal to each other. Furthermore, the strange 
C-odd component can be determined at a level of about 10%, sufficient to test for intrinsic strangeness, 
i.e. whether the C-odd component is closer in size to zero or to the C-even component. The 'instanton' 
and 'skyrmion' scenarios can thus also be distinguished at the level of several standard deviations. 

Of course, only experimental errors have been considered so far. In Ref. [^] it has been shown 
that theoretical uncertainties on first moments are dominated by the small-x extrapolation and higher- 
order corrections. The error due to the small-x extrapolation is a consequence of the limited kinematic 
coverage. This will only be reduced once beam energies higher than envisaged in this study will be 
achieved; otherwise, this uncertainty could become the dominant one and hamper an accurate determi- 
nation of first moments. On the other hand, the error due to higher-order corrections could be reduced, 
since it is essentially related to the fact that available NC data must be evolved to a common scale, and 
also errors are amplified when extracting the singlet component from NC data because of the need to 
take linear combinations of SFs. Neither of these procedures is necessary if CC data with the kinematic 
coverage considered here are available. 

4.5.2 Results for X distributions 

The best-fit SFs corresponding to the 'anomaly' refit (third column of Table |l|) are displayed as functions 



of X at the scale corresponding to the bin 4 GeV^ < < 8 GeV^, and compared with the data in Fig. 16 
Note that the SFs gi and (75 always have opposite signs because the (dominant) quark component in gi 
and g^ has the opposite sign, while the antiquark component has the same sign. For comparison, we also 
display the SFs at the initial scale of the fits, and at a high scale. The good quality of the fits is apparent 
from these plots. 

Given the poor quality of current knowledge of the shape of polarized parton distributions, it is 
difficult to envisage detailed scenarios and perform a quantitative analysis of the various shape parame- 
ters, as we did for first moments. However, it is possible to get a rough estimate of the impact of CC data 
on our knowledge of the x dependence of individual parton distributions by considering the following 
combinations of SFs, which, at leading order, are directly related to individual parton distributions: 



\ {gf- - gf-) = A. + Ac; \ {gf' + 55^^) = + Ac; (33) 

\ {aY" - aD = Ad + A.; \ {gY- + gf') = Ad + As. (34) 
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Fig. 17: The combinations of SFs of Eq. (B3|), and the corresponding parton distributions. 
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Fig. 18: The combinations of SFs of Eq. (M), and the corresponding parton distributions. 



In Figs. 17 and 18 we show, respectively, the combinations of Eqs. (33) and (P4|) for a proton target, 
together with the pseudodata for the same combinations of SFs, in the bin 4 GeV^ < <^ GeV^. In 
each figure we also display the two parton distributions, which contribute at leading order to the relevant 
combination of SFs at = 7 GeV^, as well as (as/vrjA^r at the same scale. 

Let us consider the leftmost plot in Fig. |l^. It is apparent that the expected statistical accuracy 
is very good for all data with x > 0.1. This suggests that an accurate determination of the shape of 
Au + Ac is possible. Furthermore, it is also clear that Ac (dotted curve) is extremely small with respect 
to Au (solid curve). However, we observe that the difference between the An distribution (solid) and the 
data is of the order of 15% to 20% for all x below 0.4. This difference is entirely due to NLO corrections. 
Specifically, the gluon contribution (dot-dashed curve), which spoils the leading-order identification of 
the quark parton distribution with the SF, as discussed in Sect. 4.1, Eq. (28), is small but non negligible. 
Because the various contributions to NLO corrections (in particular the gluon distribution) are affected by 
sizeable theoretical uncertainties [22], this implies that An can only be determined with an error that is 
considerably larger than the experimental one. At larger scales, the subleading corrections to coefficient 
functions are expected to be smaller and smaller, while a residual gluon contribution persists, because of 
the axial anomaly [23]. 

A similar analysis of the left plot of Fig. [I^ tells us that a determination of the shape of An is 
essentially impossible. This combination of SFs is the preferred one for a determination of the charm 
distribution, since perturbatively we expect Ac = Ac, and An is much smaller than An. Nevertheless, 
it is apparent from this figure that even in this case a determination of the charm distribution is out of 
reach. 
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A study of the down quark and antiquark distributions can be similarly performed by looking at 



Fig. |17[ The conclusion for Ad is similar, although perhaps slightly less optimistic, to that for Ati: a 
reasonable determination of its shape is possible, but with sizeable theoretical uncertainties. Likewise, 
the conclusion for As is similar to that on Ac, namely, a determination of its shape is out of reach. The 
lower plot shows that no significant information on the shape of Ad or As can be obtained from this 
analysis. An alternative handle for the determination of As and As is provided by the study of events 
with a tagged final-state charm, as discussed in the next subsection. 



4.6 Extraction of As(x) and As(x) from tagged charm 

At leading order, charm quarks are produced by CC DIS off a strange or a down quark. The combination 
of strange and down quark distributions is determined by the CKM quark-mixing matrix and can be 
written as an effective distribution 

s^(x, Q2) = \V,s\h{x, + \V,dfd{x, Q2) , (35) 

where \Vcs\'^ 0.95 and \Vcd\'^ ^ 0.05 are the squared CKM matrix elements. The same expression is 
valid for the antiquark distributions and for the polarized distributions. 

For heavy-quark production, the x and determined from the momentum of the final-state lep- 
ton cannot be directly related to the momentum fraction and virtuality of the scattered parton. Taking 
kinematical corrections into account, one finds | |3T| ] that the parton distributions are probed at momentum 
fraction ^ = x(l + mf./Q'^) and virtuahty ^u^ = + m^. 

The cross section for charm-quark production in CC DIS can be obtained from the inclusive cross 
section ([I]) by replacing the inclusive SFs with charm production ones, which read at leading order: 

fZ\x,q') = s^(e,/x2), 



W+(^ r\2\ _ ocr^ctc ,,2\ 



FZ~{x,Q^) = -2s-(e,/x2). (36) 

An alternative procedure to incorporate effects from the charm-quark mass into the cross section has 
been proposed in [^2|]. The predictions obtained in either prescription differ by less than 15%. 

The CC production of charmed hadrons in the final state off polarized targets allows a direct 
measurement of the polarized strange quark and antiquark distributions. The polarized cross-section 



difference can be obtained from the expression for inclusive CC DIS (gO|), by replacing the SFs; the 
leading-order expressions read: 



<75^/(x,Q2) = -As^(e,/i2^ 



"cJi 



gZ~{^,Q^) = As-(e,//2). (37) 
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Fig. 19: Estimated statistical errors on the polarized strange and antistrange quark distributions, compared to the GS(A) [^] 
and GRSVv |]34| parametrizations. 



The resulting leading-order cross-section asymmetry ( 1 8 ) is the ratio of polarized to unpolarized effective 
strange-quark distributions: 



(38) 



With and s^{£,,fic) being known from high-statistics measurements off unpolarized targets, a 

measurement of these asymmetries can be turned into a determination of the polarized strange quark and 
antiquark distributions. 

The statistical error on an extraction of the polarized distributions from these asymmetries is given 

byd: 



(39) 



where N^^y is the number of CC charm production events per target polarization. The estimate of sta- 
tistical errors follows the analysis detailed above. The target correction factors which account 
for the target density, the incomplete target polarization and the dilution factor (cross-section-weighted 
fraction of polarizable to unpolarizable nucleons in the target), have to be re-evaluated. Taking into ac- 
count the charm production cross sections off protons and isoscalar nucleons, we obtain F*^* = 3.3 and 
ptgt _ 2^ assume 100% charm reconstruction efficiency. It is expected that the experiments to be 
built for a j^-Factory will have efficiencies rather close to this optimal choice. 

In estimating the statistical errors, we use the same specifications as in the inclusive studies on 
polarized targets above: E"^ = 50 GeV, A'^^ = lO^*^ per target polarization, decaying in a straight section 
of 100 m length, a proton target with target density of 10 g/cm^, a detector with radius of 50 cm positioned 
at distance of 30 m from the end of the straight section. The cuts applied on the events are a minimum 
cut on the final-state muon energy of 3 GeV and minimum cut on the partonic centre-of-mass energy 
equal to the charm-quark mass. We use the same bins in x and as in the inclusive studies (Fig. ||), 
and compute the weighted mean values of ^ and /x^ for each bin. 

Figures [l9|a,b show the expected statistical errors on As'^(^, ^1) and As'^(^, /i^). Since the polar- 
ized antidown-quark distribution is not expected to be substantially larger than the polarized antistrange 
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Fig. 20: Polarized heavy quark production cross section at NLO, using the polarized GRSVv [ p4t parton distribution functions. 

one, one can determine As(^, /x^) ~ As^(^, /x^). The effective polarized strange-quark distribution does, 
however, receive a significant contribution from the polarized down-quark distribution. Approximating 
the relative error on A.d{x, Q"^) to be 10% over the full range in x and (Fig. 18), we can estimate 



the error on As(^, /i^) extracted from As'^(^, /x^)- The result is shown in Fig. [I9| c. The comparison of 
Figs. [T9| b and |T^ shows that a possible discrepancy between As(^, ^u^) and As(^, /i^) (as suggested for 
the unpolarized distributions in [^) could be detected at a level of about 20%. 

A major uncertainty in the extraction of the polarized strange-quark distribution from charm-quark 
production arises from higher-order QCD corrections, consistent with the fact, discussed at the end of 
Sect. 4.1, that the singlet quark distribution is affected by large factorization scheme ambiguities. The 
NLO contribution from the boson-gluon fusion process to heavy-quark production is proportional to the 
size of the polarized-gluon distribution, which is at present only constrained very loosely from the scale 
dependence of the inclusive polarized SFs. Figure ^ illustrates the relative magnitude of leading and 
NLO quark and gluon contributions for the GRSVv polarized parton distribution. It can be seen that 
the next-to-leading order gluon-induced subprocess amounts to a 50% correction for this distribution. 



It follows that the NLO error estimates of Figs. |17yi8| cannot be compared directly to the LO error 
estimates of Fig. which do not include the uncertainty from higher order gluonic contributions. For 
other parameterizations of polarized parton distributions the effect is smaller (since the strange quark 
distributions are in general assumed to be larger than in GRSVv), amounting typically to a gluonic 
contribution of 25%. Note also that to NLO a finite renormalization is necessary in order to relate 
quark distributions given in the MS scheme (such as GRSVv) with those of the AB-scheme shown in 



Figs. |17|,|i8[ This transformation is given for first moments in Eq. (|29|). 

As in the case of the results obtained form global fits of inclusive SFs, discussed in the previous 
subsection, the appearance of the gluon contribution at NLO poses the most significant limitation to the 
extraction of the polarized strange densities, even when using tagged-charm final states. An accurate 
knowledge of the polarized gluon density is therefore a mandatory ingredient for the fuU statistical po- 
tential of the z^-Factory to be exploited. The COMPASS experiment will provide a measurement of the 
polarized gluon distribution in the kinematic range relevant to the present studies. To which extent this 
new knowledge will improve the prospects for the extraction of the polarized strange component of the 
proton, will be known once the COMR^SS results are available. 



5 MEASUREMENTS OF as 

The value of the strong coupling as is one of the fundamental parameters of nature. There is al- 
most no limit to our need to determine it with more and more precision. The study of scaling vio- 
lations of DIS SFs and the deviation from the quark-parton model prediction of DIS sum rules, e.g. 
of the Gross-Llewellyn Smith (GLS) sum rule [35], have provided in the past, and still provide, an 
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important framework for measuring Og. Good examples are the recent NLO global analysis [^] of ex- 
isting charged-lepton DIS data, giving as(Mz) = 0.1165 ± 0.0017(stat + syst) (theor) and 
the detailed NLO studies of the 1993/94 HERA data for Ff at small x and large Q^, giving 
as{Mz) = 0.122 ± 0.004(exp) ± 0.009(theor). Theoretical errors of the latter value of as{Mz) 
are dominated by the renormalization and factorization scheme ambiguities, and by ambiguities in the 
resummation of small-x logarithms. The former might be reduced after taking into account next-to-next- 
to-leading order (NNLO) perturbative QCD effects. Indeed, NNLO combined fits to the charged-lepton 



DIS data [ ]37| ] give smaller uncertainties. The latter could be reduced thanks to recent theoretical progress 



in the small-x resummation [38]. 



A NNLO analysis of the uN DIS data of the CCFR collaboration for xF^ ||3g| gives as(M^) = 
0.118 ± 0.002(stat) ± 0.005(syst) ± 0.003(theor) [||] (more detailed studies of these data are now 
in progress pT]]). This value should be compared with the independent NLO extraction of from the 
xFs and F2 data of the same collaboration: as{Mz) = 0.1222 ± 0.0048(exp) ± 0.0040(theor) 



Other available estimates of as{Mz), including accurate NNLO results obtained from the analysis of 



LEP data and of r decays, can be found in recent extensive reviews ||43|, In this section we will 
study the potential impact of future SF measurements at the i/-Factory on the determination of as- The 
influence of the higher-twist (HT) terms, which become important in the extraction of as at relatively 
low energies, will be analysed as well. 

5.1 Determination of Og and higher-twist terms from QCD fits of the SF data 

To estimate the projected uncertainties on the determination of as at the z^-Factory, we applied the same 



procedure as we used in the estimate of the uncertainties in PDFs (see Section [3.2| ). The power corrections 
of 0(1/(5^) were also included into the generation of the fake i/-Factory data and in the fits. The target 



mass corrections (TMC) of 0{M /Q ) were taken into account following Ref. [l45[]. The additional 



dynamical twist-4 non-perturbative contributions were parameterized in the additive form: 

F2{x,Q') = F2{x,Q')^™^ + H,{x)^^, (40) 



xFsix,Q') = xFsix,Q')^™^ + Hsix)^-^, (41) 



LT TMC 

where -F2 3 results of NLO QCD calculations with TMC included, and are parametrized 

at X = 0, 0.1, ... , 0.8 and linearly interpolated between these points. Since the fitted functions depend 
on H2^3 linearly, the errors on the coefficients of H{x) at x = 0, 0.1, . . . , 0.8 do not depend on their 
central values. 

The NLO fit to the generated F2 and XF3 "data" returns a statistical error Aas{Mz) = 0.00029. 
This is much better than the statistical error on ag obtained in the global analysis of charged-lepton DIS 
data [§] and it is 16 times smaller than the one obtained in the analysis of the CCFR neutrino DIS data 



of Ref. [p%, which used a model-independent description of the HT effects []42|]. 

We verified that the extraction of the error is quite stable against changes in the PDF parametriza- 
tion. To obtain this result we repeated the analysis, generating central values of F2 and XF3 and using 



the parametrization of these SFs obtained from the NLO fit to the CCFR data given in Ref. |42]. Even 
though the functional form and the number of F2 and XF3 parameters used in the two cases are quite 
different, the difference in the errors on as is negligible with respect to the statistical accuracy of the 
individual fits. 

We also carried out a NLO fit using only the XF3 data. This results in Aas {Mz) = 0.00074, which 
is about 2 times larger than the uncertainty obtained from the fit to the combined F2 and X-F3 data. It must 
be pointed out, however, that the use of F2 in the fit introduces a strong correlation between the value of 
as and the sea and gluon densities, leading to a potential source of further systematics. In addition, the 
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Fig. 21: The errors on HT contributions accessible at the z^-Factory. 

fit to XF3 only has the advantage that inclusion of higher-order QCD corrections will be simpler, since 
the NNLO corrections to the coefficient function of the DGLAP equation are known [^]. Moreover, a 
model of the NNLO non-singlet (NS) splitting function [47, ^] was also recently constructed, using the 



results of the explicit analytical calculation of the NNLO corrections to the NS anomalous dimensions, at 



fixed number of NS Mellin moments pQ , plQ and using additional theoretical information, contained 



inRefs. [52 



I (for more details see for instance the review []54|]). For this reason the analysis of xF^ 
may give useful information on ag both at NLO and at NNLO. 

The expected precision of the i/-Factory data exceeds that of available measurements of DIS SFs 
at moderate Q^. This may allow us to improve our knowledge of the HT contributions i^2,3- The errors 
on the coefficients of the functions //2,3> which can be obtained from the analysis of z^-Factory data, are 
given in Fig. These errors were obtained from the fits described above, alongside the errors on as 
(the simultaneous estimate of the HT and 0:3 errors is very important in view of possible correlations 
between them). The errors on the HT contributions are smaller by over one order of magnitude than 
those extracted from the CCFR data [@||]. 

The data obtained at the i/-Factory could then be used for the verification of the models describing 
the HT terms. One of these models is based on the application of the infrared renormalon (IRR) technique 



57]. The advantage of this approach is that it connects the HT contributions with the leading twist 



ones. For example, in the NS approximation the IRR model of twist-4 contributions has the following 
form: 



H2,3{x)=A^I dzC2,3{z)F^J{x/z,Q) 



(42) 



Here 6*2,3(2;) are calculated in Ref. [58], and the parameter A2 introduced there can be expressed as 



(43) 



where Cp = 4/3 and /3o is the first coefficient of the QCD /3-f unction. (A similar definition was used in 
Ref. [59] in the comparison of the IRR model predictions for F2,3,L with the data.) The parameters 3 



for the IRR model can be extracted from the i/-Factory data with the errors AA| = 0.0030 GeV^ and 
AAi = 0.037 GeV^ (for comparison, the preliminary results of the combined analysis of the CCFR [ ]39[ ] 
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and the JINR-IHEP 



collaborations data are A3 



0.44 ± 0.19 GeV^ and = 0.91 ± 0.77 GeV^ 



It is worth stressing that the results of the NNLO fits to the CCFR data [40, 55], which use the 



NNLO QCD expressions for the Mellin moments of xF^ and the anomalous dimensions calculated in 



Refs. [49, 50], demonstrate the effect of shadowing of the twist-4 terms by higher-order perturbative 
QCD corrections. It should also be stressed that the decrease of the size of the fitted HT contributions at 
the NNLO level is confirmed independently by the DGLAP analysis of the combined F2 charged-lepton 
data made by the MRST collaboration [^3|], which incorporates NNLO corrections to both the coefficient 



function [ |64| ] and the model of the splitting function [)47[]. However, since Ref. [ ]63[ ] did not assign 
errors to the HT terms extracted at NLO and NNLO, we cannot decide whether the effects observed 



in Refs. [ ]40| , p5p arise from the incorporation of the NNLO corrections into DIS fits, or whether they 
demonstrate a lack of precision of the analysed data. The possible analysis of more precise data from the 
i/-Factory may allow us to clarify this point. We evaluate that the correlation coefficients between as and 
if2,3 are not so large: their maximal value is ~ —0.7 at x ~ 0.5. This allows the unambiguous separation 
of the logarithmic-like and power-like contributions to the Bjorken scaling violation. In particular, it is 
possible to hope that a clearer separation of the twist-4 effects from the perturbative QCD contributions 
may be possible at the NNLO level. The detailed study of this problem is rather intriguing. 



5.2 Determination of as from the Gross-Llewellyn Smith sum rule 

The value of as can be also determined from the GLS integral 



(44) 



At 0{al) and including the 0{1/Q'^) corrections, the GLS integral for / = 4 massless active flavours 
is equal to 



3.25 r^ii^^lV- 12.2 ^^^^MV' 



vr 



TT 



h 



(45) 



The GLS integral has been measured in a number of experiments (see for instance Ref. [ ]660 for a review). 
Its dependence was extracted by combining the CCFR data for XF3 [|^] with those from CERN and 



IHEP experiments, at several energy bins 

In many other processes the theoretical as uncertainty is dominated by the error due to the trun- 
cation of the higher-order perturbative QCD corrections. Since the theoretical expression for the GLS 
integral is known up to 3-loop as corrections, the scale- and the scheme-dependence ambiguities on the 



as extraction from the GLS sum rule can be minimized (see Ref. [ 68]). The mentioned theoretical uncer- 
tainties can survive if the Q^-region of the data used for the estimate of the integral is large and the 0(a^ 
approximation needs to be used to interpolate data from different Q'^ regions [|69|]. This problem can be 
avoided if the data are split into relatively small bins, as is expected at the z^-Factory. An additional 
contribution to the GLS sum rule comes from heavy quarks (c, s). The heavy-quark mass correction is 
known at ©(og) [[70|]. Its effect is small at energies close to the threshold, and is comparable in size with 
estimates of the massless 0(ol^ correction made with different methods [^T]]. Together with the mass- 
less contributions, the mass-dependent terms are therefore under control. In the asymptotic regime, these 



affect the threshold matching conditions [|72[], and introduce an uncertainty of about 6.5mg in the choice 
of the matching point. This uncertainty leads to an additional theoretical ambiguity of approximately 



0.002 on the value of a^{Mz) (see e.g. Ref. [gg]). 



While completing our report we learned of the new xF-^ data, obtained recently by the Hi Collaboration at HERA ||62|]. 
These data are related to rather high region ( >1500 GeV^). We therefore expect no essential improvement of our 
estimates from the inclusion of HI data into these fits. 



30 



Table 2: The statistical errors on GLS integrals at different Q bins, obtained from the different data sets (I; generated data for 
//-Factory only; II: the same data for the //-Factory combined with the CCFR data of Ref.[E^). 
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Fig. 22: The outer error bars are determined by the GLS errors of column I of Table g, while the inner ones are fixed by the 
GLS errors of column II. The points at largest Q^-bins are extracted from the CCFR data only. The curve shows the errors on 
as due to the uncertainties of higher-twist contributions. 

An important source of experimental uncertainty on the measured GLS integral is the error due 
to the extrapolation of xFs to the unmeasured high- and low-x regions. This error can be large if the 
neutrino energy is limited, as is the case for the 50 GeV //-Factory option. To estimate this error for 



J pi 

xF:i{x) = ^x''{l-xf , A= dxx''~^{l 
A Jo 



(46) 



different values of we split the xF^ data, used in the analysis of Subsection 5J^, in several bins of Q^, 
and then generate random xF^ values in each bin with the central values given by 

A- JO 
The statistical errors are given by the study of Section ^. The values of the parameters used for the 
generation were chosen as /a = 3, a = 0.7, 6 = 4. Then Eq. (^ is fitted to the generated data in each 
bin with the parameters a, b, and I3 set free. The uncertainty in the fitted value of I3, which gives 
the uncertainty on the GLS integral, accounts for the uncertainty in the extrapolation to the unmeasured 
X regions. The errors on the fitted values of are given in Table ^ One can see that at high the 
errors are quite large. Combining the i/-Factory data with the CCFR measurements of Ref. [ |39| ] one 
can obtain a significant improvement in the precision of the GLS integral determination, since the two 
sets of data have a complementary x coverage. This is shown in the second column of Table |2[, where 
the errors are significantly smaller than those obtained in Ref. [67] from the analysis of the combined 
CERN-FNAL-IHEP data with similar binning in Q^. 

The estimates of the uncertainties on as, which can be obtained from the measurements of the 



GLS integral at a //-Factory, are given in Fig. 22. The central points for as shown in this figure were 
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calculated from the solutions of the 2-loop and 3-loop renormalization group equations with the boundary 
value as{Mz) = 0.118. The error bars of the points were obtained by rescaling the statistical errors on 
/a given in column II of Table ^ by the factor of dScLS / da^ in each bin. The statistical error on as 
is smaller at 0{a'l) since \dS/das \ is larger at this order. The errors given in Fig. 22 are small enough 
to allow for the clear observation of the dependence of the QCD coupling constant measured in one 
single process and the comparison with various theoretical predictions. 

In order to estimate the as(M^) precision accessible from the GLS measurements, we fitted the 



expression of Eq. ( p5| ) to the points of Fig. 22. The analysis was made both at 0{as) and dX 0{ 



Since in our analysis the low-Q data from the i/-Factory are used, an important source of the total ag 



error is related to the error in the HT parameter h of Eq. (|45|). It should be stressed that, contrary to the 
x-dependence of {{^{x), the theoretical estimates of its first moment, related to the GLS sum rule, are 
theoretically more solid. In fact in this case we know the expression for the local operator contributing 
to the dynamical 0{1/Q'^) correction [73]. Moreover, there are several model-dependent calculations 
of the value of this operator. The first one comes from the QCD sum rules method [[74|], implemented 



in its 3 -point function realization in the work of Ref. [ J/SQ and later on in Ref. [1760. Another estimate 
of the value of the twist-4 contribution to the GLS sum rule comes from the instanton vacuum model 
[[TT]]. Within theoretical errors, the results of the model-dependent calculations of Ref s. [75, 7^ 77 1 are 
in agreement. 

To extract as{Mz) with a model-independent treatment of the HT corrections, we fitted the ex- 
pression of Eq. ( p31 ) with the parameter h set free. At ©(og), and using the analysis of the GLS data with 
the errors from column II of Table ^ we get the following results: 



At 0{as) we get instead: 



Aas(Mz) = 0.0035, A/i = 0.13 GeV^ 



Aas(Mz) = 0.0039 , Ah = 0.07 GeV^ . 



(47) 



(48) 



The cause of the only marginal improvement in accuracy when going to higher order is the faster decrease 
of ScLS with Q^, and the consequent increased correlation of ag with the HT term h/Q'^. The influence 
of the uncertainties of the high-twist correction on the precision of the Os determination at various 



is illustrated in Fig. 22, where the value Ah/Q rescaled with the factor dScis/das is given as well. 



Notice from Fig. |22| that this uncertainty weakly depends on the perturbative approximation for the GLS 
sum rule used in the analysis. One can also see that at small the errors on as due to HT uncertainties 
are increasing. Since the values for as at large have large statistical errors, the related value of 
as{Mz) is mainly determined by the uncertainties of HT corrections and does not change significantly 
from the 0{as) fit to the 0{al) one. 

If one fixes h, the statistical error on ag {Mz) in the 0(ag) fit to the data with errors from column II 
of Table ^ reduces to 0.00026. In this case, however, the uncertainty on ag due to the model dependence 
of h 0] is large (see e.g. For this reason, the error on Og obtained by considering the 



model-dependent estimates of the HT contributions to the GLS sum rule is essentially the same as the 
one defined from the existing neutrino DIS data. Therefore it seems more appropriate to analyse the 
GLS sum rule data from the i/-Factory using the fit with the model-independent definition of the HT 
contribution. 



5.3 Measurement of Fi (x) and unpolarized Bjorken sum rule 

While Fi{x,Q'^) is known theoretically to be related to F2{x) via the Callan-Gross relation and its 
calculable higher-order corrections, only recently have the experiments attempted a direct measurement 
from the data. Preliminary results on the determination of Q^) from the large y = Ehad/Ey 

behaviour of d^avN /dxdy have been obtained by the CHORUS collaboration at CERN [78] and by the 
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Fig. 23: Theoretical errors on as{Mz) from the i^-Factory measurement of the unpolarized Bjorlcen sum rule at different Q 
values. The dotted line correspond to the contributions of the TMC uncertainty; the dashed one to the uncertainty due to the 
twist-4 contribution; the solid line is the combination of both. 

CCFR-NuTeV collaboration at Fermilab [79]. However, only a few points for Fi{x, Q"^) over a limited 



range of x were extracted up to now. As discussed above, the large statistics available at the z^-Factory 
will in principle allow a complete separation of the various SF components, including in particular a 
measurement of Fi{x,Q'^). An example of the accuracy with which the individual components will 
be extracted was given in Fig. ^ This improved knowledge on Fi can be used for different purposes, 
and in particular for an independent measurement of as- This possibility is based on the study of the 



unpolarized Bjorken sum rule, which was derived in Ref. []80|]. This old, but still experimentally untested 
NS combination of neutrino DIS SFs, has the following form: 



Si 



dx 







pun 



pup 



which in terms of parton distributions can be expressed as 

»i 



n—p 



dx 







n(x, Q') - u{x, Q') + d{x, Q') - d{x, Q') 



(49) 



(50) 



Taking into account the correction of 0{as) calculated in Ref. [ ]8T| ] and twist-4 terms, the theoretical 
expression for Si~^ reads 

3 vr + Q2 • 



gn-p 



(51) 



The massless corrections to 5" ^ of 0(al) were calculated in Ref.[82], while the 0{al) contributions 
are analytically evaluated in Ref. [Q]. The heavy-quark mass correction to S^^^ is known from the 
calculations of Ref. [ |70| ] and is comparable with the results of existing 0{ag) estimates Jti]]. 

The HT term in Eq. ( |5T| ) is analogous to that of the GLS sum rule. The value of h bj is proportional 
to the matrix element of a local twist-4 operator, hBj = — (8/9)((C')) [f73|], with: 



I a 



and Of, = uGf^vlvlbU - dG^^vlul^d 



(52) 



where G^^y = (e^,ya/3/2)G^^(A"/2). The application of the 3-point function QCD sum rules results in 



0.15 lb 0.07 GeV [75], where we take for the theoretical error the 



the following estimate: {{O)) 
conservative estimate of 50%. 

We constructed the 0{as) evolution of (see Eqs. (p^,(|50|) using the set of parton distri- 
butions of Ref. [|8|], and taking into account the twist-4 contribution, and estimated the theoretical error 
on as{Mz) extracted from this expression. The behaviour of the error Aq!s(M^) that can be obtained 
from the measurements of the 5f " integral at different is shown in Fig. This error is defined as 

-^(Q2) dasiQ^y-' 



dS' 



das(g2) das{Mz) 



(53) 
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where AS'f " is the error due to the different sources of theoretical uncertainties. One of them comes 
from the error in the contribution of the TMC, which is related to the uncertainty in the existing sets of 
PDFs. We calculated this error using the uncertainties of PDFs from Ref. ^ and convinced ourselves 
that it does not exceed 0.002 (see Fig. |23|). At O (as) dS^~^ /das = 2/ (Stt), and the theoretical error due 
to the HT uncertainty is {3tt/2) x O.OT/Q^. For the reference scales = 4 GeV^ and =iq GeV^, 
we obtain A^-^as(M^) = 0.012 and A^'^as{Mz) = 0.007, respectively. To keep a balance between 
the various sources of uncertainties the corresponding statistical errors must not exceed these values. 



From the theoretical point of view, the uncertainties of Eq. ( p3| ) are consistent with those of 



asiMz) extractions from the GLS sum rule value at = 3 GeV^ [68], namely A^^ as(Mz) = 0.003. 



The latter one is a bit smaller because of the differences between the perturbative expressions for the GLS 
sum rule and the unpolarized Bjorken sum rule, and between the estimates of the twist-4 contributions 
to these sum rules. It is therefore rather important to check the results for the high-twist contributions 



to Si ^ obtained in Ref. [75]. Within the framework of the instanton model this question is now under 
studyQ. 

The experimental determination of the Bjorken unpolarized sum rule, which will be possible at 
the z^-Factory, can therefore be considered as an additional source for the determination of as, provided 
the twist-4 contributions are known with more precision. 

Precise z^-Factory data on Fi{x, Q"^) will also allow the measurement of the x-dependence of the 



HT contribution to Fi, which in the IRR model of Ref. [ jSSQ is predicted to coincide with the shape 
of the twist-4 contributions to xF-^. In view of the considerable interest given to the analysis of the 
contributions of HT terms to different quantities, it is quite desirable to study this prediction in detail, 
using the experimental data for Fi. 

To conclude this section we note that the measurement of the unpolarized Bjorken sum rule re- 
quires the extraction of F^^^x, Q"^) from DIS on a hydrogen target and of Q^) from DIS on a 
deuterium target. The latter process necessitates the analysis of nuclear corrections, especially in the 
small-x region. A detailed study of these problems, as well as the discussion of other experimental 
alternatives, will be discussed in the next section. 

6 NUCLEAR EFFECTS IN DIS AT THE z^-Factory 

There are two general motivations to study nuclear effects in DIS experiments at the i/-Factory. First, 
nuclear physics of parton distributions is of interest in itself, and the comparison of heavy-target data 
with hydrogen and light-nuclei data (e.g. deuterium) may give us new insights into the structure of 
multi-quark systems. On the other side, an accurate knowledge of nuclear effects is necessary in order to 
extract the SFs of a physical proton and neutron from nuclear data. This applies in the first place to the 
neutron, since available neutron targets are mainly nuclei. Theoretical studies of nuclear effects, among 
other possible applications, could also help in choosing the most appropriate neutron target. 

DIS from different nuclear targets has been studied with electromagnetic fi/e probes at CERN, 



SLAG, and FNAL (for a recent review and references, see [|84|, |85|]). It was observed that heavy-target 
SFs differ substantially from those of light nuclei in a wide kinematical region of x and Q^. Figure U 
presents a compilation of data on the so-called EMC ratio (a traditional measure of the magnitude of 
nuclear effects in DIS), F^/F^, where F^ and F^ are the SFs per nucleon of a nucleus with mass 
number A and of deuterium, respectively. One passes through several distinct regions with characteristic 
nuclear effects when going from small to large x. At x < 0.1 one observes a systematic reduction of 



the nuclear SFs, the so-called nuclear shadowing. This is illustrated in the right-hand panel of Fig. g4 
showing the EMC ratios on a logarithmic scale. A small enhancement appears there at 0.1 < x < 0.3, 
followed by a dip at 0.3 < x < 0.8, which is usually refen^ed to as the 'EMC effect', and finally an 
enhancement, which is associated with nuclear Fermi motion. 



''C. Weiss, private communication. 
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Fig. 24: The x dependence of heavy-target/deuteron SFs ratio as measured in DIS of muons and electrons off various nuclear 
targets and averaged over Q^. The left panel shows data from CERN and SLAC for different nuclear targets. The right panel 
focuses on the region of small x and illustrates nuclear shadowing effect. Data points are from NMC ||86|], BCDMS [p7||, SLAC 
and and FNAL (E665 collaboration |^9|) with only statistical errors shown. Nuclear targets are specified in the brackets 
on the plot legends. 



Experimental information about nuclear effects in other DIS observables, such as the ratio of 
longitudinal to transverse cross sections or the spin SF gi, is available but scarse. We note also that DIS 
off nuclear targets is characterized by additional (new) SFs, which do not appear in DIS off an isolated 
nucleon. As an example we refer to the tensor SF bi, which is specific for spin-1 targets and appears in 
DIS on deuterium (for a review and references see for instance [|^]). 

6.1 Nuclear shadowing 

Before we turn to the discussion of the DIS regime, it is useful to discuss the low-Q^ region away from 
scaling. Here the behaviour of neutrino cross sections (SFs) is quite different from that of charged lep- 
tons. For the latter it is well known that the longitudinal SF F^, as well as F2, vanish at low Q^, because 
of electromagnetic current conservation. It was shown long ago by Adler that, at low Q^, CC neutrino 
interactions are dominated by the axial current, and neutrino cross sections can be expressed through 
PCAC in terms of pion cross sections []90|]. In contrast to charged-lepton scattering, F£ is finite, domi- 
nated by a pion pole for at the pion mass scale, and drives the neutrino cross section in this region. 
Using the Adler relation. Bell predicted nuclear-shadowing effects for neutrino scattering similar to what 
is observed in pion-nucleus interactions [^I]]. Going to larger brings a finite contribution from vector 
and axial-vector meson states, which have been discussed in terms of an extension of the vector meson 



dominance model to vector and axial- vector currents [ p2| , p3| ]. Charged-current neutrino interactions 
with nuclear targets were studied in bubble-chamber experiments [^], where nuclear shadowing was 
observed at low Q^. 

Most of the attempts to understand nuclear shadowing are based on the space-time picture of DIS 
at small x in the target rest frame, where DIS is viewed as the process of interaction of the partonic (or 
hadronic) component of the exchanged 7* or W* with the target. At small x the typical propagation 
length of those states exceeds the average distance between bound nucleons, and coherent effects in 
the propagation of partons through the nuclear medium are important. Nuclear shadowing is usually 



explained by multiple-scattering effects from bound nucleons []85|]. 

Nuclear shadowing in F2 was calculated for both low and high regimes in terms of two differ- 



ent models in [95] in an attempt to match muon (NMC) and neutrino (CCFR), data [[39|] on F2 at small 
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X. It was found that nuclear shadowing in F2 is similar (though slightly smaller in magnitude) to that 
observed in muon-induced reactions. 

At large in the scaling regime both charged-lepton and neutrino-induced reactions are described 
by universal parton distributions. Some observations on nuclear modifications of different combinations 
of parton distributions can be made from existing charged-lepton DIS and Drell-Yan data. Phenomeno- 
logical constraints on the behaviour of nuclear sea and valence quarks at small x were discussed in 
[100, 101 , 102|]. Explicit evaluations of nuclear effects in singlet and non-singlet combinations of parton 



distributions were performed in [ |103| ], where nuclear shadowing for F2 and xF^ was studied in terms of 



the non-perturbative parton model of [ |104[ ], which was extended and applied to nuclear targets in [ ]105| ]. 
It was found that, while the shadowing effect in neutrino F2 is similar to the corresponding effect in 
charged-lepton DIS, the nuclear shadowing for xF^ is enhanced with respect to that for F2 in the region 



of small X < 0.01 (see Fig. g5|). It was argued in QlOSO that the underlying reason for the enhancement 
of nuclear shadowing for xF^ is its negative C-parity. In the small-x region, xF^ is determined by the 
difference of effective quark and antiquark cross sections, and it is known from the multiple-scattering 
theory that the double-scattering correction to the difference of the cross sections is up to a factor of 2 
larger than the corresponding correction to their sum. 



Nuclear shadowing effect for F2 and F3 




0.0001 



0.001 



Fig. 25: The ratios of a heavy target to the free nucleon SFs R2 = F'{^ /F-i' and Rs — / F^^ calculated for the ^^Fe nucleus 



in the region of small s at = 10 GeV^ niO: 



We note in this respect that a similar enhancement of nuclear shadowing was predicted for the spin 
SF gi (see discussion in [|^]), which involves the differences of quark and antiquark distributions with 
helicities parallel and antiparallel with respect to the helicity of the target. 



6.2 Nuclear effects at large x 

The physics mechanisms that generate characteristic nuclear effects at large x are quite different from 
those that govern nuclear shadowing at small x. At large x the typical DIS time scale in the laboratory 
reference frame is small with respect to an average distance between bound nucleons. This allows us to 
assume that nuclear DIS is dominated by incoherent scattering from bound nucleons. It was found long 



''This disagreement has recently been resolved by CCFR/NuTeV [ Pq | who employed, among other things, a proper treatment 
of the charm mass threshold effects [0, ^. The ratio of the new F2 values measured in and fi scattering is now in agreement 
with the NLO predictions, which use the massive charm production scheme implemented in the MRST parton distributions 
setl|9il. 
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ago that major nuclear effects here are due to nuclear binding [ ]105| , |106i |107| ], which leads to a depletion 



of nuclear SFs at x ~ 0.5, and to the Fermi motion [ |108[ ], which is responsible for the enhancement at 
X > 0.7. These effects explain the bulk of the observed behaviour of nuclear SFs at a; > 0.2, though 
detailed understanding of this region is far from complete and further studies of the reaction mechanism 
are required. 

It is quite important to separate nuclear effects in a QCD analysis of neutrino data. As an example 
we refer to the recent analysis of higher-twist terms in the CCFR data on the xF^ SF [ ]109| ]. The theoreti- 
cal and experimental situation becomes less clear when going to the region of x close to 1. Here we enter 
into the resonance region, and the notion of twist expansion becomes less well defined. Nuclear effects 
in SFs are essential in this region, because the nucleon inelastic SFs must vanish as x —>^ 1. The impact 



of nuclear effects on HT terms in xF^ was studied in [ ]109| ], where it was shown that the consideration 
of nuclear effects at large x somewhat decreases the magnitude of dynamical HT terms extracted from 
CCFR data on iron target. However, the results of calculations of nuclear corrections are sensitive to the 
details of the nuclear structure input, e.g. the behaviour of the nuclear spectral function in the region 
of high excitation energy of the residual nuclear system, which are not known yet. In order to minimize 
uncertainties associated with nuclear effects, it is desirable to use light nuclear targets in the experiments 
at the z/-Factory. 

We note also that nuclear SFs can extend beyond x = 1, the kinematical limit for scattering from 
a free nucleon. Events with x > 1 have indeed been observed for F2 in fi DIS from a carbon target by 



the BCDMS collaboration 1 110] and recently by the CCFR collaboration in neutrino DIS from an iron 
target [HI]. It is rather interesting to search for similar events at the z/-Factory, where the statistics will 
be much higher. 

The region of large x will be explored in more detail in electron scattering experiments at Jefferson 
Lab. However the experiments at the zv-Factory at large x with different nuclear targets are challenging 
and could significantly contribute to the field by providing a direct measurement of the EMC effect for 
different parton combinations, such as the C-even F2 and C-odd xF^. 

6.3 Nuclear effects in DIS sum rules 



As discussed in Section |5.2| , the GLS sum rule is a convenient tool to extract as from neutrino data. 
Neutrino data are usually collected on heavy nuclear targets; therefore, it is of importance to separate 
contributions to the GLS integral associated with nuclear effects. We denote Sq-^^ = S'qlq + 6Sgls> 
where Sqlq and Sqi^^ are the GLS integral for the nucleus of A nucleons and the isoscalar nucleon 
respectively, and 6Sgls accumulates corrections due to nuclear effects. There is a number of effects 
that can contribute to 5Sgls- Contributions due to nuclear binding, Fermi motion and off-shell effects 



were discussed in [ ]112| ]. It was found that these effects cancel out in the leading twist, and a tiny 
correction appears as a higher twist. For example, 6SqI^^ = —1.2 x 10~^GeV^/(5^ and SSqi^^ = 
-1.9 X lO-^GeVVQ^ for the iron and deuterium nuclei respectively. These quantities are more than an 



order of magnitude smaller than the corresponding QCD power correction estimated in [J75|]. However, 
as was discussed in [ |103| ], nuclear shadowing gives a finite and rather large negative correction already in 
the leading twist. In particular, a ~ 4% renormalization of the GLS sum rule due to nuclear shadowing 
was found for the ^Spe nucleus, ^ii-S'gls/S'^ls = -0-035, at = 10 GeV^. 

At this point we must mention that it is usually believed that the GLS sum rule is not renormalized 
by nuclear effects in the leading in l/Q"^ order, since to this order the GLS integral counts the baryon 
number of the target. However, a negative sign of the shadowing correction is also a generic feature 
of multiple scattering theory. It is therefore challenging to look for a dynamical mechanism that would 
compensate a negative nuclear shadowing correction in the GLS sum rule. Certainly more work is needed 
to clarify the status of the GLS sum rule, as well as other DIS sum rules, in nuclear targets. 

The Bjorken sum rule for the non-singlet combination of the neutrino SFs Fi was discussed in Sec- 
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tion 53 as an alternative tool to extract Og- Neutron data are necessary in order to measure this sum rule, 
and it is clear that, since there exists no free neutron target, nuclear data have to be used. Combined hy- 
drogen and deuterium data are usually used as a source of information about the neutron. The integrated 
difference between the hydrogen and the deuteron SF could then be used to extract Sf~^ = S^~^ + 6 Si, 
where the last term incorporates nuclear effects. Using the method described in |112|, it can be shown 
that the corrections due to nucleai" binding and Fermi motion cancel out in the leading twist, similar to the 
GLS sum rule, and the corresponding power l/Q^ correction is small, of the same order of magnitude 
as in the GLS sum rule. 

Note also that the deuteron nucleus, in spite of a weak binding, might not be a perfect source of 
information about neutron SFs, especially at very large or very small x. It was emphasized recently that 
DIS experiments with mirror ^He and nuclei, which form an isotopic doublet, could give information 
about the neutron SF practically free from contamination from nuclear effects [113]. We then suggest 
that the direct measurement of the difference Si ('^He) — Si (^H) could be a better source for the Bjorken 
sum rule than the corresponding hydrogen-deuterium difference. The calculation of corrections to the 
Bjorken sum rule, as well as to the more fundumental Adler sum rule, due to meson-exchange currents 
in nuclei and nuclear shadowing, is under wayQ 

In summary we note that because of a larger number of observables, which can be accessed with 
u and u beams, DIS studies at the z^-Factory can give unique information about the structure of hadrons 
and nuclei, information that is not accessible with fi/e machines. In particular, a unique opportunity 
offered by the z^-Factory is a direct measurement of nuclear sea and valence quarks distributions in a 
wide kinematical region. 



7 ELECTROWEAK STUDIES AT THE i/ FACTORY 

Experiments with z^e^ and uN have played a fundamental role in establishing the SM. Before the start 
of LEP, the best determinations of the electroweak mixing angle came from neutrino experiments, and 
even now the results from neutrino DIS at NuTeV play an important role in global analyses. 

The characteristics of a z^-Factory are such that very precise tests of the SM and its extensions 
may be possible from neutrino-electron scattering and neutrino induced DIS. As a first approximation, 
the information available from these experiments can be parametrized in terms of the uncertainty in the 
determination of the sine of the Weinberg angle, s^, = sin^ 9w- In fact, radiative corrections enter the 
two processes in different ways and one should look at these experiments as complementary measure- 
ments, akin to the present determinations of sin^ and Myy at e+e^ and hadron colliders. Moreover, 
a very precise low-energy determination of s^/ would test a different variety of new-physics scenarios 
than usual colliders. 



7.1 ue scattering 

ue' scattering provides a particularly clean probe of the electroweak coupling. There are several pro- 
cesses which contribute at a z^-Factory: 

(NC) i/^e"^f^e", P^e~^i?fj,e~ , (54) 
(NC + CC) i/ee"^fee" , Pee"^i>ee", , . . . , (55) 

(CC) i/^e"^z^eAi~ • (56) 

Events originated by or in a fi^ beam without a muon in the final state cannot be disentangled at a 
z^-Factory and must be considered together. 

^S.A. Kulagin, work in progress. During the course of writing this report we learned about the paper [ |ll4| ], where the 
nuclear shadowing corrections to the Gottfried sum rule for ^He - mirror nuclei were discussed. 
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In the case of ve ve processes, numerical values for the total cross sections are given by 



a{ye-^ ve') = 1.72 x 10~*^cm"^ x Ey[Ge\] x 



2 I 1 2 

9l + ^9r 



(57) 



where g^ ji = s^r or g^ ji = ±1/2 + s^r according to the process. 

Despite the very small cross section, the use of a 2 ton dedicated, fully active target-detector made 
of liquid CH4 [^] or of a 20 ton liquid argon [Q] time projection chamber should provide around 10^ ue 
events/year with a /i+ beam and about half of it with a fi^ beam. In this subsection we use default beam 
specifications, a 20 cm detector radius, s'^, = 0.2314, and / Ldt = 8.6 x 10^^ cm~^. 

The signal is a forward electron track with no hadronic activity and energy above a threshold Emin- 
The transverse momentum of the outcoming electron is very small, pt ~ ^/m^E^. In the configuration 
considered here the transverse momentum due to the intrinsic divergence of the beam is even smaller. The 
main source of background is quasi-elastic V(.N scattering, which can also produce a forward electron, 
but is characterized by pt ~ yJuiNEy and can be distinguished if the pt resolution of the detector is 
good. The signal-to-background ratio is expected to be better than 5 at a 50 GeV -Factory [Q], leading 
to a minor dilution of the sensitivities considered below. The beam has the advantage that quasi- 
elastic VeN scattering produces positrons instead of electrons, so that this source of background can be 
removed. 
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Fig. 26: (a) statistical uncertainty (in units 10"'') in the extraction of sin^ 6w from ue scattering as a function of 
electron energy, (b) impact of luminosity measurement at the level of 10"'^, 5 10"'', 1 10~* on the same sin^ 9w 
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The statistical sensitivity to Sy^^ is shown in Fig. ^a as a function of Emin for the /i^ and /i+ 
beams. When only integrated cross sections are considered, the /x^" beam allows for a superior 
resolution, very close to lxl0~^, while the ji^ beam is much less sensitive because of a cancellation 
among different terms. Assuming that a measurement of the outgoing electron energy E^. is possible, one 
can study the E^. dependence and achieve a better resolution in the /x^ case, (^s^ 2 x 10^"^. This is 
close to having an optimal observable. For small Emin, the sensitivity is a mild function of the threshold 
energy. 

A significant problem in ue scattering is the normalization of the cross sections. As can be seen 



from Fig. |2gb, in order to preserve the statistical sensitivity one would need a luminosity determination 
at the level of 10^^. The goal at a i^-Factory would be to reach a precision on the flux of 10^'^; this, 
however, seems inadequate. A realistic possibility is to normalize the ve rate to up.^ , namely to the 
muon regeneration process (see Eqs. (|55|),(^6|)), which occurs for Ey > m^/2me = 10.8 GeV with 
cross sections given by 
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Again, the background from v^N scattering should be considered, but it can be drastically reduced by 
a cut on the transverse momentum [^]. This mechanism could provide a normalization accuracy of 
3.6 X 10^^, but would be available only for the fi'^ beam. On the other hand, it could serve to calibrate 
different luminosity measurements. 

We have also investigated the case in which the muon beams are polarized and concluded that 
muon polarization would not add significantly to the electroweak measurement. At best, polarization 
asymmetries could help overcome the normalization problem. With 70% polarized muon beams the 
statistical error is 5s^/ 3 x 10~^ for the |U~ beam, and worse for the /i+ beam. On the other hand, 
uncertainties in the polarization of the fi beam do affect the final precision. In order to maintain the final 
precision at the level of 6syy ^ 10^^, it is necessary to know the polarization of the muons to 0.05% or 
better. For this reason precision electroweak measurements require a storage ring design that minimizes 
these uncertainties. 

From the theoretical point of view, all one-loop QED and electroweak corrections are known 



[115] and could be easily implemented. The uncertainty from higher-order electroweak corrections can 



be certainly brought below Ss'u, = 1 x 10 [1 16]. QED corrections to the e , ^ spectra are relatively 



important and may need consideration of some higher order effect. A major theoretical uncertainty is 
likely to come from the hadronic contribution to the amplitude, which must be calculated at the 

relevant « 10^^[GeV]£'jy, a region where perturbation theory cannot be applied. The problem is 
analogous to — but should not be confused with — the one of the hadronic contribution to the running 
of the electromagnetic coupling, Aa^adr^ which enters most electroweak tests. Aahadr is needed to 
relate the fine-structure constant a{q^ = 0) to a(M|), which is relevant at the pole and for the 
electroweak corrections to muon decay as well as for low-energy NCs. The latter, however, have an 
additional sensitivity to hadronic loops in the 7-Z° mixing. This contribution can be calculated from 
e+e~ data using dispersion relations, SU(3) flavour symmetry, and perturbative QCD. The most recent 
estimate [117] leads to 6s'^, = 5 x 10^^. In view of recent progress [ ]118[ ], this can probably be reduced 



by a factor of 2 or more. However, the use of SU(3) flavour symmetry implies a sizeable ambiguity, 
which cannot be resolved by better data only. For what concerns the uncertainty in Aahadr, it should 
not be considered a limiting factor for this experiment, as it affects in the same way most electroweak 
observables. Moreover, there has been and there will be progress in its determination, and it is likely to 
play a lesser role than the uncertainty from the 'y-Z^ mixing in ue scattering. 

In conclusion, a total uncertainty on s'^, of about 2 x 10^^ is probably achievable at a z^-Factory, in 
the case where high-performance detectors are available, the polarization of the muons can be controlled 
very precisely, and progress is achieved in estimating the hadronic effects. Higher precision would 
require a substantial increase in luminosity as well as major theoretical improvements, mainly in hadronic 
physics. 

7.2 sin^ from DIS 

Current electroweak analyses of vN DIS (at NuTeV) are based on the Paschos-Wolfenstein (PW) rela- 
tion, 

which is designed to isolate the u, d valence quark contributions that cancel out in the ratio and is there- 
fore quite insensitive to the hadron structure. On the basis of about 1.3 (0.3) million (z/^) events. 



NuTeV has found [119] 



s^{OS) = 0.2253 ± 0.0019(stat) ± O.OOlO(syst) ± 0.00025(Mi) + 0.0005 In ^ . (60) 

As can be seen, if the on-shell definition s|(/(OS) = 1 — M^/M| is used, the indirect dependence 
of the result on the top and Higgs masses through radiative corrections to Eq. (^Tj) is relatively small. 
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Observable 


Stat, error 


PDF 


R^'~ 


0.4 


~ 12 




0.5 


15 


R^"' - 0.8R^'^ 


2.2 


~ 2 


P 


4.9 


~ 4 



Table 3: Uncertainties on the s'^- determination from different observables in u DIS in units 10 (see text). 



This is why the NuTeV result is often presented as an Mw determination with 6Mw ~ 130 MeV. This 
interpretation can be highly misleading when the experimental accuracy reaches 0(10"'*). 

At a z^-Factory, a striking improvement of the statistical error can be expected, as well as the 
elimination of the main systematic problems of NuTeV (uncontrolled u^. beam contamination and NC 
event identification). Whether the final s^r sensitivity may reach the level of 10^^ or dMw ~ 5 MeV 
depends, however, on many different factors. 

If the detector can identify primary in the final state, the CC events originated by electron 
neutrinos can be distinguished from NC events. One can then consider ratios of NC/CC total cross 
sections [||] 

Using our default beam specifications and Ei^pt > 3 GeV, Eh > 1 GeV, = sin^ 9w = 0.225, we 
find i?^ Rf^^ = 0.36. Also the s^^ sensitivity of the two beams is very similar, dR/dsy^ k. —0.5. 
Because of the 10^ CC and 3 x 10^ NC events available for each beam, the statistical error on the 
determination from Eq. (^) is negligible, as shown in Table ^ However, one should take into account 
that i?^ have a higher sensitivity to hadronic physics than the PW relation. As an illustration, we show 
in Table |3| the sensitivity of the determinations from R^ on present PDFs. The values are obtained 
by comparing results for i?^* for which different sets of the CTEQ5 PDFs [^] (sets m, d, and hq) have 
been employed. They are the result of a LO analysis based on present day information. Clearly, future 
improvements in our knowledge of the hadron structure at the z^-Factory and elsewhere, as well as a full 
NLO implementation of QCD radiative corrections will lower this uncertainty, but it seems unlikely that 
they will bring it down to the level of the statistical error. 

On the other hand, the PW relation is formally recovered in the combination 

{l+gr)R>'' -{r + g)R^'\ (62) 

where g = {Ey^) / {E^^) = 0.857 takes into account the different mean energy of the neutrino and 
antineutrino beams, r = crcc(j^)/o'cc('^)> and lepton universality has been used. Unfortunately, because 
of the numerical closeness of R^ and i?^^ , this combination is not an efficient probe. Moreover the s^/ 
sensitivities of R^ are also very similar and the improvement in the PDF sensitivity is paid for by the 
lower statistical sensitivity. One possibility is to fit a parameter a in i?^ — aR^^^ in order to minimize 
the overall uncertainty. For instance, a 0.8 leads to (^^^^^(PDF) w 5s|(/(stat) 0.0002. 

Of course, one can construct other ratios of cross sections, also combining the cross sections of the 
two beams in a way similar to Eq. (|^); however, they tend to be less sensitive to and/or to depend 
even more than i?^* on the hadronic structure. A possible exception could be the quantity 



(63) 



which is modelled on the PW relation and has indeed a lower sensitivity on hadronic physics (see Ta- 
ble H). However, this observable is penalized by a lower resolution and requires the knowledge of 
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the ratio acc {,IJ^^)I (^cc (/" ) — and therefore of the relative normalization of the two beams — to better 
than 10-"^. 

These high-precision measurements cannot be carried out without a high-performance tracking 
target. The detector and beam requirements based on Eqs. ( |6T1 ) have been discussed in The main ones 
are high electron/muon detection efficiency, good primary lepton identification, and precise momentum 
and energy measurements. A rough estimate of the corresponding uncertainties is (5s^ 2 x 10^^. 
The ability to have both /i+ and fi^ beams would be an important asset, as it would allow the use of 
Eq. (H). 

From a theoretical point of view, it will be necessary to implement full NLO QCD and QED 
corrections, to incorporate possibly complete two-loop electroweak corrections, to assess the importance 
of higher twists, and to take into account charm production corrections. This is to be contrasted with the 
obsolete and incomplete implementation of QCD and electroweak corrections of the NuTeV analysis. 
On the other hand, almost all theoretical tools are in principle already available, including the dominant 
two-loop weak corrections [ |116| ]. Charm production, in particular, is likely to be a significant source of 
uncertainty; it accounts for about 3% of the total CC cross section and therefore even a 1% error on the 
overall charm yield would induce a 10~^ uncertainty on i?^ or ^s^/ 2 x 10~^. An extrapolation 
from the CCFR measurements suggests (5s^ ~ 0.0003 [||]. 

If the outcoming electrons cannot be resolved and CC events initiated by electron neutrinos are 
considered as NC events — a situation analogous to present-day vN experiments — the only observable 
for the /x^ beam is 



This observable and the analogous one for the /U"*" beam have a marginally lower statistical sensitivity 
and a much higher dependence on the SFs, (5s^ > 0.02, than ii^ . 

In conclusion, efficient electron identification seems to be crucial for a very precise determination 
of Sy^r at the i/-Factory . The statistical uncertainties are likely to be negligible in comparison with 
theoretical and experimental systematic errors, which in turn are difficult to estimate at the moment. 
Before a realistic assessment of the potential of u DIS at a i/-Factory for precision electroweak tests is 
possible, several systematics will have to be understood; most importantly the precision of the i/-Factory 
measurement of the unpolarized SFs, primary lepton identification, isospin-violating effects and charm 
production. 

7.3 New physics through radiative corrections 

The very precise determination of the electroweak mixing angle from z/e and uN scattering at a z/-Factory 
would test the SM at a level competitive with LEP, SLD and Tevatron measurements. Experiments at 
the front end of a muon storage ring could therefore severely constrain many extensions of the SM, 
which potentially affect the SM predictions for neutrino scattering through virtual contributions. If the 
new physics is characterized by a high mass scale and it affects dominantly the two-point functions, its 
contributions can be parametrized in a model independent way by the 5, T, U amplitudes introduced in 



[120]. Similar strategies have been presented by a number of authors (see for instance [121]). Many 



models of physics beyond the SM can be studied at least approximately in this simple way. 

We use Mj{ ~ 100 GeV and s'^r — sin^ 9]y[AIz)y^ — 0.2315 as reference SM values and set 
g = {Ey^) / {Ey^) = 0.857. As far as ve scattering is concerned, in the case of the /x^ beam we can 
consider the ratio of electron- and muon-production total cross sections 

r- = -) { . (65) 

Virtual effects parametrized by S, T, and U shift the measured value of r with respect to the SM predic- 
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tion ro by 

r- = (1 + 0.0055 5 + 0.0004 T) . (66) 

This observable is therefore very insensitive to the isospin breaking parameter T, similarly to the weak 
charge in atomic parity violation experiments with cesium. It could allow for a very accurate determi- 
nation of the isospin conserving parameter S, within ±0.1; this is roughly the accuracy of the present 



global fit (see Ref. [ 122], p. 105). A similar observable can be constructed for the ^+ beam, provided the 
luminosity determination is based on a CC process analogous to muon production. In this case the shift 
is 

r+ = r J (1 + 0.0111 S - 0.0100 T) , (67) 
and implies a very different constraint in the (5, T) plane. Finally, the measured value of the Paschos- 



Wolfenstein relation (59) is affected by the oblique parameters according to 

i?- = i?Q (1 - 0.0136 5" + 0.0253 T) (68) 

The different ue and uN measurements would therefore provide complementary constraints on new 
physics. 



8 STUDIES WITH HEAVY QUARKS 

As has been pointed out in the past, and as shown in a previous section, the i^-Factory can be seen as 
a z^DIS charm factory. A comprehensive review of the possible physics goals of such a facility can be 
found in Ref. [||]. Here we shall explore in quantitative detail some interesting example, drawing from 
the expertise available within the Working Group. In particular we shall focus on two aspects of charm 
production in neutrino interactions: low-multiplicity processes such as the diffractive Ds and the quasi- 
elastic charmed-baryon production. We shall show that these processes allow a clear identification of the 
charmed hadron, and therefore a very good estimate of absolute decay branching ratios. Together with 
low systematic errors, these measurements can provide, for instance, a precise measurement of fos- 

In the following we consider nuclear emulsions both as neutrino target and tracking device. It is 
worth stressing that the use of nuclear emulsions is limited by the overlapping of interactions. A density 
of interactions of about 20 per cm'^ is reasonable. On the other hand, 10'' interactions are needed for 
the measurements we will discuss in the following. Such a statistics could be obtained by running the 
machine at low luminosity and taking data for a few years (months) whether the experiment is located 
far from (0(1 km)) or close to (0(100 m)) the neutrino source. 

Our simulations include realistic estimates of the experimental efficiencies and systematics. It is 
important to point out, however, that the methods discussed below can be used by electronic experiments 
with a very good vertex detector as well, although with different efficiencies and backgrounds. 



8.1 Direct evaluation of the A+ branching ratios 

8.1.1 Model-independent extraction of BR{A'^ pK~Tr^) 

So far, only model-dependent extractions of A+ branching ratios have been obtained, see [122]. A 
method, based on the neutrino quasi-elastic charm production, for a model-independent determination of 
most of the A+ branching ratios has been proposed in Ref. [123]. So far, two different methods to extract 
A+ branching ratios have been used [ ]122| ]. As discussed in Ref. [122] they rely on different theoretical 
assumptions on B physics, namely the B branching ratios to A^ and give results that are not in quite a 
good agreement. Therefore, a model-independent determination of A+ branching ratios would provide 
a better theoretical understanding of the baryonic 6-decays. For a detailed discussion of this method, see 
Ref. [pi. 
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cr(10-4o cm2)\ Model 


F.R. [121 


S.L. [1251 


A.K.K. [126|, 127, 128 j 


A.G.Y.O. [129]] 


K. [Oojl 




0.2(0.030%) 


9.0(1.3%) 


8.0(1.2%) 


1.0(0.15%) 


3.0(0.45%) 




0.6(0.089%) 


16.0(2.4%) 


10.0(1.5%) 


0.6(0.089%) 






1.0(0.15%) 


23.0(3.4%) 


41.0(6.1%) 


3.0(0.45%) 


5.0(0.74%) 




0.1(0.015%) 


5.0(0.74%) 




0.6(0.089%) 


1.5(0.22%) 




0.3(0.045%) 


8.0(1.2%) 




0.3(0.045%) 




Total 


2.2(0.33%) 


61.0(9.0%) 


59.0(8.9%) 


5.5(0.82%) 


9.5(1.41%) 



Table 4: Predicted quasi-elastic charm-production cross section assuming a neutrino energy of 10 GeV. In brackets the quasi- 
elastic charm rate with respect to deep-inelastic interactions is given. 



Experiment 


(10-« cm2) 


(10-40 cm2) 


(10-40 cm2) 


(10-40 cm2) 


(10-40 cm2) 


[ 


131 


1 


(2.3t?-]) 










[ 


132 


1 


(2.3 ±2.0) 


(4.5 ±4.0) 








[ 


133 


1 






(0.9tr7) 






[ 


134 


1 






(3.7t^-i) 






[ 


135 


1 








(38.3 ±23.1) 


< 8.8 


Average 


(2.3 ± 1.5) 


(4.5 ±4.0) 


(1.1 ±0.8) 


(38.3 ±23.1) 


< 8.8 



Table 5: Summary of the experimental measurements of quasi-elastic charm-production cross sections. 



8.1.2 Present knowledge, theoretical and experimental, of neutrino quasi-elastic charm production 

The simplest exclusive charm-production reaction is the quasi-elastic process where a d valence quark 
is changed into a c quark, thus transforming the target nucleon into a charmed baryon. Explicitly, the 
quasi-elastic reactions are 



/i"A+(2285) 



;U"S+(2455) ,v^n ^-S*+(2520) , 



i/^p ^ /i-S++(2455) ,v^p 



/u-i;*++(2520) . 



(69) 
(70) 



For a detailed theoretical review of neutrino quasi-elastic charm production, see Refs.[ ]124| ]-[ ]130| ]. 

The predicted cross sections, assuming a neutrino energy of 10 GeV and according to different 
authors, are shown in Table 0. As we can see, these predictions can even differ by one order of magnitude. 
If we express the total quasi-elastic charm production rate with respect to deep-inelastic CC interactions, 
it ranges from 0.33% to 9.0% (see Table g). From Fig. 3 of Refs. [ [I25| ] and [ pOl ] we note that the total 
cross section in both models is almost flat for neutrino energies above 8 GeV. 

The statistics of neutrino quasi-elastic charm events collected by bubble chamber and emulsion 
experiments is rather poor. The measured cross sections, obtained by a reanalysis that uses the latest 
results on A+ branching ratios 1 122 ], are shown in Table |[ Despite the large statistical error, these 



measurements are clearly inconsistent with the predictions of Refs. [ ]125| ]-[ ]128[ ], while the agreement is 
fair for the ones in Refs. [ 124, 129, 13(]| ]. The average value of the cross sections predicted by Refs. | |124 , 
129| , |130| ] has been used to get a rough estimate of the expected number of events, given in Table g 



8.1.3 Description of the method 

The method consists in identifying the Ac by means of the peculiar topology of the quasi-elastic reaction. 



As shown in Fig. g7|, the only charged particles produced in this process are the muon and the short-lived 
particle (Ac) except for the reaction Fig. p7|c) where an additional vr is produced. Therefore there is only 
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Reaction 


cr(10-39 cm2) 


n(%) 


Expected events 




0.14 


0.14 


14000 




0.06 


0.06 


6000 




0.3 


0.3 


30000 




0.07 


0.07 


7000 




0.03 


0.03 


3000 


All 


0.6 


0.6 


60000 



Table 6: Quasi-elastic charm-production cross section and its contribution to the total CC neutrino cross section. The last 
column shows the expected number of events, assuming a starting sample of 10^ CC neutrino-induced events. 



a small contamination of these events from deep-inelastic charm production where a charmed hadron is 
produced (faking a Ac) and the topology is quasi-elastic-like. 

In this way, since we only use topological information, the Ac identification is model-independent. 
The relative contamination of and Df from deep-inelastic events, ejake^ can be dealt with as a 
relative systematic error on the branching ratios. We assume the relative systematic error to be e fake + 
3cr/afce where a fake is the error on Efake ■ 

The normalization to determine the A+ absolute branching ratios is simply given by the number of 



events with a vertex topology consistent with Fig. |27|. No model-dependent information is used to define 
the normalization. The little knowledge we have about the quasi-elastic charm-production cross section, 
which is model-dependent unless measured, plays a role only in the evaluation of the deep-inelastic 
contamination, namely the systematic error. In fact, in the evaluation of the contamination, the ratio TZ 
between charm quasi-elastic and standard deep-inelastic production appears. It is worth noting that, even 
if the ratio IZ had an uncertainty of 500%^, the relative systematic error on the branching ratios would 
be ~ 7.2% (see [ll23|]). 






a) 



b) 



c) 



Fig. 27: Topology of the quasi-elastic neutrino-induced charm events in the case of the reactions a) i^^n ^ /i~A^, b) u^n — + 
/i~Sj(E*^) and c) — > /i^Ej^(Sc^+), with subsequent Ec decay into Ac. The particles inside the box represent the 
A J decay products. 



8.1.4 Measurement accuracy 

The expected accuracy on the determination of the A+ branching ratios as a function of the relative error 
on TZ, the quasi-elastic charm production cross section relative to the deep-inelastic one, is shown in 
Table ^ To compute the expected number of events in each decay channel we use the central values 
(shown in Table ^ together with their errors) given by the Particle Data Group [122]. From this table we 



Nevertheless, this is not the case. From Table H we see that IZ is measured with an accuracy better than 100%. 
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Channel 


PDG BR [122] 


ABR 


ABR 


ABR 






= 10%) 


= 100%) 


= 500%) 


A+ pK^TT+ 


(5.0±1.3)% 


(±0.09±0.04)% 


(±0.09±0.09)% 


(±0.09±0.4)% 


A+ ^ Afi+Vf, 


(2.0±0.7)% 


(±0.06±0.01)% 


(±0.06±0.04)% 


(±0.06±0.1)% 


A+ ke+Ve 


(2.1±0.7)% 


(±0.06±0.01)% 


(±0.06±0.04)% 


(±0.06±0.1)% 



Table 7: Statistical and systematic accuracy achievable in the determination of the absolute branching ratios, assuming a 
collected statistics of 10^ CC events, as a function of the relative error on TZ. The central values are taken from Ref. 1 122]. 



can see that, even assuming a very large (unrealistic) systematic error, the achievable accuracy on the 



branching ratios makes the discrimination between the methods discussed in Ref. [122] still possible 



8.2 Direct evaluation of branching ratios and io^ measurement 

The experimental knowledge on leptonic Dg decays is very Uttle. Currently, the branching ratios for 



Ds lu decays have been estimated by the PDG | [T22| ] to be BR{Ds fii^) = (4.6 ± 1.9) x lO'^ 
and BR{Ds —>■ tv) = (7 ± 4) x 10^^. These large uncertainties translate into a large uncertainty on the 
extraction of the decay constant fo^ ■ 

A method that would allow at the z^-Factory the extraction of most of the Ds branching ratios. 



and consequently of //j^, by means of purely leptonic decays, has been proposed in Ref. [ |136| ]; the 
expected accuracy would be better than 5%. Once fo^ will be measured with such an accuracy, one will 
feel more confident about extrapolating to the decay constants in the B system, fs and fs^, which are 



crucial quantities for the quantitative understanding of B^gy 
(Vts) from them. 



B^^s^ oscillations and the extraction of Vtd 



8.2.1 Topology of neutrino-induced diffractive charm events and background 

In Z)i*^[] diffractive production, only a muon is produced at the interaction point (primary vertex), besides 
the charmed meson. Therefore, these events are characterized by a peculiar topology: two charged tracks 
at the primary vertex, one of them being a short-lived particle. 

Neutrino-induced quasi-elastic charm events are characterized by the topologies shown in Fig. 27. 
Therefore, they are similar to the diffractive ones, but with a cross section twice as large. Since antineu- 
trinos cannot induce quasi-elastic charm production, while diffractive production is the same for both 
V and v, in the following we will consider only u beams. We will then make the assumption that all 
the events with the above topology are due io ' diffractive production. In this case we will wrongly 
classify some of the charmed hadrons produced in deep-inelastic interactions. 

The charm production in v interactions and the event topology have been studied by using the 



HERWIG event generator [ |137p , as well as an event generator based on JETSET [ ]138| ] and LEPTO [ [13^ 
assuming the energy dependence of the charm fractions reported in [14C]^. The average charm fractions, 
convoluted with the neutrino spectrum, are Fgo = 61%, Fj^- = 26%, F^- = 7.3% and 



a; 



5.7%. 



The signal kinematics has been studied by using an event generator developed within the CHORUS 



Collaboration [141|. 



The contamination to the diffractive sample, which comes from deep-inelastic events, can be writ- 



ten as 



(71) 



'in the following, -Oi*' means either Dg or D*. The same notation is also used for D mesons. 
*We assumed that charm fractions are equal for both v and v. 
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fake/ £ fake 


£fake{%) 


eDi'){%) 


15% 


23% 


0.037±0.009 


3.3 ±0.8 


30% 


34% 


0.04±0.01 


3.3 ±0.8 


50% 


53% 


0.04±0.02 


3.3 ±0.8 


100% 


101% 


0.04±0.04 


3.3 ±0.8 



Table 8: The relative and absolute error on e/afce as a function of the relative error on TZ. In the fourth column the systematic 
uncertainty due to the Z)'*' contamination is reported. 

where 

^ a{9,N ^ /x+pW- iV) 



We take the charm production in u interaction to be 3% [ |142| ]. By using the charm fractions 
discussed above, {F^- + Fj^-) = 31.7%. The factor ffake = (6.0 ± 0.1)% is the fraction of deep- 
inelastic charmed events faking the diffractive topology; Skm = (0.4 ± 0.2)% gives the efficiency of 
kinematical cuts as explained in Ref. [ ]136| ]. Finally we get e/afce — 0.04%. 



The currently measured TZ has an error of about 15% (see Ref. [136]), which affects the relative 



error on ejake- As described in Ref. [ |136| ], a kinematical analysis allows a good Dg detection efficiency 
(~ 73%) with a small background. The expected contamination and its error are given in Table |8[ 

Another possible source of irreducible background is the diffractive production of D^*'>~ (see 
Ref. [136]). It is suppressed by a factor |Fcs/K:d|^ ~ 20 and only 30% of the D*^ decay into a charged 



charmed meson {BR{D*- D') = 0.323 ± 0.006 [p2|]). On the other hand, all the diffractively 
produced are background: we assume that are half of the diffractive sample ^ Finally, we get 
ED* = (3.3 ±0.8)%. 

From the numbers given above, it turns out that the little knowledge we have about the diffractive 
charm production cross-section plays a role only in the evaluation of the deep-inelastic contamination, 
namely a term of the systematic error. Even if the ratio IZ had an uncertainty of 100%]^ the relative 
systematic error on the branching ratios would be ~ 0.04% (see Table |8|). We want to stress that, since 
the contribution is dominant, the overall systematic uncertainty does not depend at all on the IZ 
accuracy. Writing the relative systematic error from deep-inelastic events as ejake + '^'^fake^ where a fake 
is the error on efake-> the overall relative systematic uncertainty is 



£sys = ^{e fake + ^'y fake? + £%.) = (3.3 ± 0.8)% , (73) 

which is dominated by the ££)(,) term. 
8.2.2 Description of the method 

An almost pure sample of Z)^ from diffractive events, with a small contamination of and pro- 
duced in deep-inelastic events and in diffractive production, can be built by using diffractive Di*-* 
production from antineutrinos. The normalization to determine the Dg absolute branching ratios is given 
by the number of events with a vertex topology consistent with one ^ plus a short-lived particle. No 
model-dependent information is used to define the normalization. 

It is worth noting that the contamination of and A^ events does not affect the Dg ^ t channel. 
Indeed, such events would present a unique topology with two subsequent kinks. An event with a double 



kink has been observed recently in CHORUS (see Ref. [ |141| ]). 



This assumption has been drive n by the NuTev results: oiv^iN ^ DsN) — (1.4 ± 0.3) fb/nucleon, a{Ufi,N 



li'DlN) = (1.6 ± 0.4) fb/nucleon [ |143| ]. 

'"Nevertheless, this is not the case. We recall that the Big Bubble Chamber Neutrino Collaboration 1 144 ] and NuTeV 1 14; 
combined analysis gives an accuracy of about 15% for 1Z. 
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Channel 


PDG BR [12211 


New method 




(4.6 ± 1.9) X 10^3 


(±0.55±0.15) X lO-''' 


Ds TV 


(7 ± 4)% 


(±0.17±0.23)% 


Ds (plv 


(2.0 ±0.5)% 


(±0.08±0.07)% 



Table 9: Statistical and systematic accuracy achievable in the determination of the Ds absolute branching ratios, assuming a 



8.2.3 Measurement accuracy at a neutrino factory 

At present there are no experiments with both an adequate spatial resolution to fully exploit the diffractive 
topology and a sufficient antineutrino-induced CC event statistics. Therefore, the method proposed in 
this paper could only be exploited with the above-mentioned detector exposed at a i/-Factory. 

Let us assume that 10^ v CC events are collected into an emulsion target and that the detection 
efficiency is about 73% for the Dg decays. By assuming a vertex location efficiency of about 50%[3 and 
assuming a v diffractive production rate of 6.2 x 10^'^/CC event, we expect to detect a number of Dg 
equal to Nd^ = lO'^ x 6.2 x lO^^ x 0.73 x 0.5 ~ 2.3 x 10^. 

The expected accuracy on the determination of the Ds branching ratios is shown in Table ^ for a 
few channels, together with the current status. To compute the expected number of events in each decay 
channel we have used the central values (shown in Table ^together with their errors) given by the Particle 
Data Group JT22| ]. 



As discussed in Ref. [ |145| ] the leptonic branching ratios are proportional to the decay constant. 
Therefore, by using the measured branching ratios given in Table ^ /^i^ can be extracted. If we collect 
10^ u^jL CC interactions we get fn^ = 288 ± 4(stat) it 5(syst) MeV, where the central value is taken 



from Ref. [122] 



8.3 Theoretical estimates for z^-induced exclusive production 

The observation of exclusive Ds production at the i^-Factory, through the process v^+N^ fi'^ +N+D~ , 
is also of interest for the study of the production mechanism within QCD and opens a new possibility to 



study the nucleon structure. There exists a QCD factorization theorem [146, 147], which states that the 
amplitude for hard exclusive meson-production processes can be written as a convolution of a skewed 
parton distribution (SPD), a distribution amplitude, and a hard part. This theorem has recently been 



applied to the investigation of electroproduction of single light mesons [ |148| ] and meson pairs [|14? 
Motivated by the possible implications for the z^-Factory, the formalism has been extended to CC-induced 
processes]^. We summarize here the main results of this work, whose validity is limited to values of 
large compared to —t and to the squared masses of the involved particles. The differential cross section 
of the leptoproduction process is given by 



(74) 



where I is the neutrino momentum and T is the amplitude for the subprocess 

W£\q) + N{p)^D-{q') + N{p') . (75) 
At leading order, T is obtained from the sum of three diagrams involving a gluon SPD (Fig. 28 a and 



diagrams obtained by an interchange of the order of the vertices) and two diagrams with a contribution 
of the (polarized and unpolarized) strange quark SPD (Fig. |28| b plus one diagram with a changed order 
of vertices). The relative Feyimian diagrams are convoluted with the gluon and s quark SPD and with 

"This efficiency accounts for the electronic detector reconstruction and the automatic location of the event vertex inside the 
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Fig. 28: Two of the five contributing diagrams. 

the distribution amplitude (z) for the meson (z being the fraction of the D~ momentum carried 
by the strange quark). 



For a numerical estimate of the cross section we adopt the distribution amplitude from [15C] 



<^ (z) = NdJz{1- z)exp 



m 



2^2 



(76) 



where uj = 1.38 GeV was obtained in [ ]151[ ] as the best fit for the D meson, and N^^ has to be chosen to 
satisfy the sum rule 

[\z^ {z)=fD^ {11) 

Jo 

For a comparison, we shall also use the asymptotic form of the Dg distribution amplitude: 



cl>„-(z) = 6/z5,^(l-z). 



(78) 



In both cases we use fn^ = 270 MeV, namely the average of the results obtained so far in lattice 



calculations [152]. The gluon and quark SPD's are parameterized as in [149], combining Radyushkin's 
model [147, 153 ] with the parameterizations of the usual (forward) parton distributions of Ref. [ ]154| ]. 

Figure ^ shows the results obtained for the differential cross sections da/d^Bj and dcj/dQ^ 



where t = {p — p')'^ has been integrated over the interval tn 



assuming for simplicity a fixed neutrino energy Ep = 34 GeV. For the plot of the xbj -dependence 
has been integrated from 12 GeV^ to the upper bound given by the constraint y < 1, with y := 
P ■ q/p ■ J- = Q"^/ i'^^BiP • 0- The plot of dcr/dQ^ is based on the xbj -integrated cross section, integrated 
between x^j = 0.18 and xbj = 0.75 and taking into account the same kinematical constraint. The 
solid lines correspond to the form of <I>^- given in Eq. (77), while the dashed lines correspond to the 



< -t <2 GeV" 



asymptotic form of Eq. (78). The dotted lines are obtained by setting to zero the strange quark SPD, 
proving the dominance of the gluon contribution, and suggesting that this process is a potential probe for 
the measurement of the gluon SPD. 

Integrating over all variables Q^, x^j, and t over the same kinematical region gives a value for the 
total cross section of a = 5.6 x 10"^^ GeV~^ = 2.2 x 10~^ pb. Uncertainties of this rough estimate 
result from the current limited knowledge of the SPDs, by the 20% uncertainty on fn^ , and by the lack 
of knowledge about the exact form of the meson-distribution amplitude, as illustrated in Fig. 29. It 



is worth noting, however, that experiments of the kind discussed here can provide much more precise 
information on the Dg decay constant independently of the exact cross section. As discussed in the 
previous subsection the relatively high production rate of Dg mesons allows us to determine by 
measuring the Dg branching ratios and its total width [ ]136| ]. In spite of the small cross section, the huge 
integrated luminosities of / Ldt > 10^ pb^^ available at the i/-Factory would lead to samples of the 
order of 10^ events. 

The analysis shows that the large rates, and the extended range available, will allow a more 
accurate determination of the gluon SPD and help to better test the theory of these exclusive processes. 



emulsions. 

''B. Lehmann-Dronke and A. Schafer, in preparation. 
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Fig. 29: The differential cross section for exclusive D~ production as a function of xbj or respectively. The dotted lines 
show the contribution stemming from the gluon SPD. The results obtained for the asymptotic form of the distribution amplitude 
<J>^- are plotted with dashed lines. 

The available statistics may even allow a study of higher-twist corrections to the leading behaviour, 
probing the limits of the factorization theorem. 

9 A POLARIZATION IN NEUTRINO DIS 



Many experiments have reported [ ]155| ]-[ ]158| ] the observation of longitudinal polarization of A baryons 
produced in neutrino DIS on an unpolarized target. Measurements in the current fragmentation region 
(CFR) give information on distribution and fragmentation functions, while measurements in the target 
fragmentation region (TFR) allow to access fracture functions [ ]159| ]. In both cases the measurement of 
A polarization provides a sensitive way for studying perturbative and non-perturbative spin phenomena. 

9.1 A polarization in the current fragmentation region 

The polarization of spin 1/2 baryons inclusively produced in polarized DIS in the CFR may be useful 
to obtain new information on polarized distribution and fragmentation functions. A lot of theoretical 



attention has been dedicated to the self -revealing polarization of A's and other hyperons [160]-[171] 



Most papers, with the exception of Refs. [164], [167] and [17C], do not consider weak interaction con 



tributions, since there is no available experimental information. 

The NOMAD collaboration has published some interesting results [158] on the A polarization in 
CC interactions; more data, sensitive to 7/Z interference effects, might be available from high-energy 
NC processes at HERA; more complete information is however only expected from experiments at the 
z^-Factory. 

These experiments are a unique source of new data, owing to the natural neutrino polarization and 
to the selected couplings of W^'s to pure helicity states. We will be able to study in detail processes like 

vp^rk^ + X vp^tA^ + X (CC) 
up^v + X vp^D +X (NC) 

where the proton p may or may not be polarized, depending on the experimental setup, whereas neutrinos 
are obviously always polarized {\y = —1/2, Ap = +1/2). Notice that though the formulae we present 
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here hold for proton targets, they can be easily modified into analogous expressions valid for polarized 
and unpolarized neutrons. We consider CC and NC processes separately. The explicit expressions for the 
polarization of the final A's in terms of elementary dynamics, quark distribution and fragmentation func- 
tions show how these experiments can provide precious information on distribution and fragmentation 
functions which are still far from being well known. 

9.1.1 Charged current neutrino processes, vp^ih^ X 

For the neutrino-initiated processes, the longitudinal polarizations P[u/] and P[y/] for any spin 1/2 
baryon B (A's and A's for instance) produced in CC DIS processes are defined as 

Pw/iB) = ^ ^^up^e-B^x (^9) 

and 

^^up^£+B+X _ ^^up^l+B^X 
P\-^A = ^^9p^£+B+X + ^^up^i+B.X ' 

where B± denotes a baryon B with helicity it. 

In the most general case, when also the proton p is polarized - and we denote by a superscript S 
its spin state - we obtain at leading twist in the QCD factorization theorem: 

piS) ^ j:.M)-^da'^^"^DB/u. - rfgy^-AZ)^/,-.] ^^^^ 

Y.t,j[[Uir.'da^_ 'DB/d, + {dj)\'d(j^^ Db/uA 

(S) 

where (g)^ stands for the number density (distribution function) of quarks q with helicity ± inside a 
proton with spin S, whereas q± alone will refer, as usual, to a proton with helicity -i-. For the flavours we 
use the notation Ui = u,c and dj = d, s. The polarized fragmentation functions are defined in terms of 
fixed-heUcity fragmentation functions as 

^Db/q = DB+/q+ - DB-/q+ = ^B-fq- " ^B+fq- , (§3) 

and da'^Z^^^ stands for the da/dy cross section for the elementary interaction udj^iui, where ly and 
dj having negative helicities. 

The above polarizations depend on the usual DIS variables x, z and y; apart from the de- 
pendence of distribution and fragmentation functions due to the QCD evolution, there is a kind of fac- 
torization in the dependence on the three variables, in that the distribution functions depend on x, the 
fragmentation functions on z and the SM dynamics on y. If convenient, and according to experimental 
setups, numerator and denominator of Eqs. ( ^l] ) and (82) can be integrated over some variables. 



and 



Performing the sum over flavours in the numerators and denominators, neglecting c-quark contri- 



butions and inserting the elementary dynamics expressions, Eqs. ( pl| ) and (|J) give, for longitudinally 
polarized (it helicity) protons 



and 



p(±).„ , [d±+R s±] ADb/u -{l-yfu^ [ADB/d + R ^Db/s] 
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where R = sin^ 6^/ cos^ Oc ^ 0.056. 

In the simpler case in which the proton is unpolarized, we replace q± by q/2, and Eqs. (84) and 
(p^) become respectively 



[^'^1 ^ ' ' ^' ^ [d + R s] Ds/u + {l-yru [D^/, + R D^/,] ' ^''^ 

p(0) , „ , + i^g] Ai^B/^ - (1 - yf u[ADs/d + RADb/s] 

[^,^,y,^)- [ J ^ ^ ^^^^ + (1 - y)2 + R D^/,] ' 

The formulae given above hold for any baryon and antibaryon with spin 1/2; further simplifications are 
possible when a A baryon (and, more in general, a baryon rather than an antibaryon) is produced: in 
this case, in the kinematical regions characterized by large values of x and z one can neglect terms that 
contain both q distributions (in a proton) and q fragmentations (into a A) as they are both small. Then 
one simply has: 

( A; - <i ( A; ^ - , (88) 

d(±)^A \ d(0) /A \ ^DA/d + R ADj^/s 

P,-IJA:z) ~ P,-iAA;z) ~ —, (89) 

and the polarizations, up to QCD evolution effects, become functions of the variable z only, since any 
other term apart from the fragmentation functions cancels out. 



[d + Rs] ADb/u -{l-y?u [ADs/d + R ^Db/,] 



and 



Equations (88) and ( [89| ) relate the values of the longitudinal polarization P(A) to a quantity with 
a clear physical meaning, i.e. the ratio AD lyj^/ Dj^, this happens with weak CC interactions - while 
it cannot happen in purely electromagnetic DIS | |171 ] - due to the selection of the quark helicity and 



flavour in the coupling with neutrinos. A measurement of P(A) offers new direct information on the 
fragmentation process. 

Weak CCs couple to pure helicity states, and do not transfer transverse polarization from quarks to 
final baryons; therefore transversely polarized protons do not add any information. However, one might 
still have - in analogy to what happens in unpolarized A^-A^ interactions - final A's with transverse 
(with respect to the production plane) polarization. This can only originate in the fragmentation process 



of an unpolarized quark; recently, new polarizing fragmentation functions have been introduced | |172[ ] to 
describe such an effect: 

A^D^,/^{z, k^) = Dj,y^{z, k^) - Dj,i/^{z, fcx) (90) 

where k± is the transverse momentum of the A with respect to the fragmenting quark momentum. A 
measurement of transverse A polarization would give a direct measurement of such a new function: 

;j(A; z, k^) ~ — '- ■ (91) 

9.1.2 Neutral current neutrino processes, up^vA^X 

In analogy to the previous paragraph, the longitudinal polarizations of the produced baryon B is given 
by 

Plu,u] [Bj — ^^^p_,^B+X _|_ (j^fjup^uB-X ' ^^^^ 

and 

da^P~^^^+-^ — da^P^^^--^ 
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For the numerator and denominator of P[iy^^j{B) and P[p^pj{B) separately, one obtains, for a 
generic spin state S of the proton (C = sin^ Ow/^)' 



AD 



B/uj 



+ 
+ 

+ 



AD 



B/dj 



(5) 



AD 



B/uj 



(5) 



AD 



B/dj 



(94) 



and 



Df 



+ (1 - y)Hc' + {djf^ (1 - 2C)2 



D 



+ 

+ 

+ 
+ 
+ 



{S) 



Bjdi 

D 



B/uj 



{d,)f{l-yY{l-2CY + {d^)'yAC 



16^2- (1-4 C)^ 
4C2 - {d^f^ (1 - yf (1 -2Cf 

{u,)f{l-ACf-{ujf^\l-yfma 
{d,)f{l-2Cf-{d^f^ {l-y)HC^ 
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B/uj 



AD 
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AD 



B/uj 
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(95) 



(96) 



^gl] (^) = E . { 16 + (1 - yf (1 - 4 C)^ 

+ 4C2 + {d^f^ (1 - y)2 (1 -2Cf 

+ (1 - 4C)2 + {u,f^ (1 - y)2 16 

+ K)f (1 - 2C)2 + {d,f^ (1 - y)2 4C2 



B/uj 



D 



B/dj 



D 



B/uj 



D 



B/d, 



}■ 



(97) 



In the case of A (or any baryon, rather than antibaryon) production, a simple expression for its 
longitudinal polarization P can be obtained by neglecting the antiquark contributions and the terms 
proportional to sin^ Q-w- For longitudinally polarized protons in this approximation we have 



{K)^(l-8C)AZ?A/„. + (d,)^(l-4C)Ai?A/,.} 

' (A) ^ — 

E,- {(uj)^ (1 - 8C)I?A/„^, + (d,)^ (1 - 4C)Z)am} 

whereas for unpolarized protons, where q± is replaced by g/2, we obtain 

^(0) , , ^ (1 - 8 C) Ai^AM + (1 - 4 C) Ai^A/., } 



(98) 



E, {^^,- (1 - 8 C)Z?A/„, + dj (1-4 C)i?A/d, } 



(99) 



9.1.3 Present knowledge on A fragmentation functions and numerical results 

The study of A polarization gives direct access to new fragmentation functions of quarks into A; it is thus 
worth looking at the present knowledge of these functions, both unpolarized and polarized. 
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Unpolarized A fragmentation functions are determined by fitting e+e~^ A + A + X experimental 
data, which are sensitive only to singlet combinations such as {D/^/g + Df^jg). It is therefore impossible 
to separate the fragmentation functions relative to A's from those for A's in a model-independent way. 
Furthermore, flavour separation is not possible without appropriate initial assumptions. For example, 
one can assume SU (3) flavour symmetry, as in Ref. [|163|]; this leads to the simple approximation 



Dk/u = D/y/d = Da/s = ^A/s = ^x/d = ^A/s ■ (100) 

Other parameterizations rely on different assumptions: for instance, SU{3) flavour and SU (6) spin- 
flavour symmetry breakings lead in Ref. [168] to -Da/s > ^K/u = ^A/d' in Ref. 1 173[ ] to Da/s ^ 



Da/u = D^/d- 

Other examples could be given, but at this stage much more stringent data are greatly needed for 
the determination of A fragmentation functions. 

For polarized A fragmentation functions the situation is also problematic. In fact, polarized A 
fragmentation functions are obtained by fitting A polarization at LEP, that is only sensitive to non-singlet 
combinations such as AD^/g ~ ^-^A/g = ^^A/g- ^^^^ ^^^^ have direct information on the va- 
lence contributions to the polarized A fragmentation function, and the sea contributions are generated 
through evolution, for each flavour, starting from a given initial scale. It is therefore possible to de- 
termine polarized fragmentation functions for A and A separately. Unfortunately, flavour separation is 
very difficult in the polarized case, because experimental data currently available cannot even discern 
between remarkably different input models adopted for the valence contributions: for example, the ratios 
Cq{z) = ADx/q{z)/ Dx/q{z) range from values such as Cg = I, Cu = Cd = in the non-relativistic 



quark model to Cs = 0.6, Cu = Cd = -0.2 in Ref. [174], to Cs = Cu = Cd = z" in one of the 



scenarios of Ref. [ |163[ ]. The parameterizations of sea and gluon fragmentations have to be guessed by 
mere assumptions. 

The most consistent and model-independent fragmentation functions that can be derived from 
e~^e^ data are (-Da/ij + ^A/q) ^"^d' separately, AD/^^g and AD^^g. These do not allow a computation 
of P(A) and P(A). The lack of knowledge of unpolarized fragmentation functions for separate A and A 
could be partiaUy overcome by introducing new measurable quantities: 

where T = da{A)/da{A). These two quantities can be computed using the singlet unpolarized frag- 
mentation functions, and they are simply related to -P(A) and -P(A) through the factor T, which can be 
measured; notice that one always has P* < P. In Fig. ^ we show, as an example, the expected values of 
P*(A) and P*(A) obtained for typical kinematical values of the planned i/-Factory experiments, adopt- 
ing two different sets of fragmentation functions from Ref. [ ]163[ ]. Note that, in the chosen kinematical 
region, one expects T ^ 1 at large z so that P*(A) ~ ^'(A), while P{A.) could be sensibly larger than 
P*(A) and in fact comparable in size with -P(A). 

Let us finally emphasize the impact that neutrino semi-inclusive DIS (SIDIS) data would have on 
our knowledge of both nucleon distribution functions and A fragmentation functions. Once the A and 
A polarizations in neutrino and antineutrino SIDIS off a polarized target will be measured, a carefully 



combined study of equations like (|4|), (|85|), (|88|), (|89|), (|98|), (|99|) would not only allow the extraction 



of new information about polarized and unpolarized A and A fragmentation functions, but also a cross 
check of our knowledge on nucleon polarized distribution functions. This combined program would offer 
an invaluable insight in the internal spin structure of hadrons in terms of their elementary constituents. 
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Fig. 30: P*(A) (left) and P*(A) (right) for the process >/i X as a function of z, as predicted by using scenarios 2 



and 3 for the polarized A, A fragmentation functions of Ref . 1 1 63 ] . Results for unpolarized (solid) and longitudinally polarized 
(A = +: dashed; A = — : dot-dashed) proton target are shown. Kinematical conditions typical for the i^-Factory have been 
considered: = 30 GeV, 0.01 < a; < 0.7, 1 GeV^ < < 100 GeV^ W^^ > 4 GeV^ 

9.2 Models for A polarization in the target fragmentation region 

To describe the A polarization in the TFR one has to model the fracture function, namely the probability 
of finding a parton q in the target nucleon and a final hadron with given momentum and polarization. 
This problem was addressed in only few models, which we discuss below. 



9.2.1 SU (6) Quark— diquark model 

Longitudinal polarization of A's produced in DIS was first considered in [175]. After kicking out a 
left-handed u or d quark from an unpolarized nucleon, we have the following relative probabilities for 
the polarization states of the remnant diquark: 
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(103) 



where, for example, {uu)iq denotes a uu-diquark with total spin S = 1 and spin projection Sz = 0. 
It is assumed that during recombination with the unpolarized s quarks the diquark does not changes 
its polarization. Since in the naive quark model (NQM) the polarization of A's is equal to the s-quark 
polarization, the A's that are directly produced will be unpolarized. However, the final state A's may 
also be produced indirectly via electromagnetic decay of S° or strong decay of S*'s. In both cases, the 
non-strange diquark changes its spin from 1 to 0, while the strange quark retains its polarization fIT6? 
Using the SU{Q) wave functions of octet and decuplet baryons, we obtain for the A polarization: 



pun 



^A 



-0.55 
-0.05, 



(104) 
(105) 
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yielding, for an isoscalar target, 



pu{p+n) ^ pPip+n) ^ _Q_3Q^ 



(106) 



9.2.2 Meson-cloud model 

Some non-perturbative features of the nucleon structure, such as the deviation from the QCD-parton- 
model-inspired Gottfried sum rule, can be explained in the framework of the meson-cloud model. The 
pion-cloud model provides a natural explanation of the isospin symmetry breaking in the unpolarized 
proton sea. In the case of polarized DIS, the scattering on the lowest A K and S K component of the 
nucleon wave function provides a possible mechanism leading to a violation of the Ellis-Jaffe sum rule. 
The polarization of A's produced in the TFR in this model has been considered in [ |176| ]. It appears that 
A polarization is almost 100% anticorrelated with the target polarization. It is thus expected to be for 
an unpolarized target. 



9.2.3 Polarized-intrinsic-strangeness model 



The polarized-intrinsic-strangeness model (for a review, see [ 177]) qualitatively reproduces experimental 
features of the 4> production in pN annihilation. The model is based on the following major observations. 
First, the fact that the masses of pseudoscalar mesons are small with respect to the typical hadronic scale 
can be attributed to the existence of effective strong attraction in the J^^ = channel. Secondly, 
from phenomenological analyses of the quark condensates in the framework of QCD sum rules, it is 
known that the vacuum density of strange-antistrange quark pairs is comparable to the density of u and 
d quarks. It is natural to assume that the polarized constituent quark can contain an ss pair with the 
vacuum quantum numbers corresponding to a ^Pq state. Hence, in the polarized nucleon, the spin of s 
will be antiparallel to the valence quark spin, Sz{s) = —1/2 for Sz{qv) = 1/2. In Ref. [162] there was 
considered the case (in the following referred to as A) of an angular momentum projection of the ss 
pair Lz{ss) = +1 (Sz{ss) = —1). In this case any s quark in the target fragment should have negative 



longitudinal polarization, so that the longitudinal A polarization should also be negative, see Fig. 3 1 




U , s 



Fig. 3 1 : Dominant diagram for A production in thie target fragmentation region due to tire scattering on a valence u-quark. Each 
small arrow represents the longitudinal polarization of the corresponding particle. 

In the quark-parton model of deep inelastic v or iy scattering, the net longitudinal polarization Pg 
of the remnant s quark is given by 
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where Nq (Nq) is the total number of events in which a quark (antiquark) is struck, and Cs q is the spin- 
correlation coefficient. The antiquarks contribute with a negative sign because their CC weak interactions 
are right-handed. The final A polarization Pa = Dp Ps, where Dp is a. dilution factor that describes the 
spin transfer during hadronization. 

Let us consider also the scenario B, where both projections Lz{ss) = +1 and Lz{ss) = of 
the ss-pair angular momentum contribute with equal probabilities. Lz{ss) = means neglecting the 
transverse motion of the ss pair [ |178| ]. 

The correlation of the remnant s-quark polarization with that of any other struck sea quark {qsea / 
s) depends on whether they come from the same parent constituent quark. If they do, which might be the 
dominant case, then a strong spin-correlation is expected (case a). Otherwise, the correlation should be 
reduced (case b). Then for the spin-correlation coefficients we have 
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The results for the remnant s-quark polarization and the predictions of NQM are presented in the Ta- 



ble 10 together with the measured A polarization. The data on longitudinal A polarization from the 



Table 10: A polarization in the TFR of (anti)neutrino SIDIS. Model predictions in boldface lie within ±1(t of the experimental 
data. 



Experiment (Reaction) 


Data 


Aa 


Ab 


Ba 


Bb 


NQM 


WA21 [ 


155 


] (l^fM - P) 


-0.29 ±0.18 


-0.51 


-0.75 


-0.22 


-0.25 


-0.55 


WA21 [ 


155 


] (i^M - P") 


-0.38 ±0.18 


-0.85 


-0.92 


-0.30 


-0.31 


0.03 


WA59 


156] (i)^ - Ne) 


-0.63 ±0.13 


-0.82 


-0.91 


-0.29 


-0.30 


-0.30 


E632 [ 


151] (u^-Ne) 


-0.43 ± 0.20 


-0.70 


-0.84 


-0.27 


-0.27 


-0.30 


NOMAD 


[158] {u^ - C) 


-0.21 ± 0.04 


-0.59 


-0.80 


-0.24 


-0.27 


-0.30 


NOMAD [ 


158 


] (^^M - 'V') 


-0.29 ± 0.06 


-0.54 


-0.77 


-0.23 


-0.26 


-0.55 


NOMAD [ 


158 


] {Jy^ - "n") 


-0.16 ±0.05 


-0.61 


-0.81 


-0.25 


-0.27 


0.03 



NOMAD experiment [ |158p have the best statistical accuracy, and we will base the conclusions mainly 
on comparisons with these data. As one can see from Table |l^: 

• the predictions of NQM are not so different from the NOMAD data on the isoscalar (mainly 
carbon) target, but contradict the NOMAD p and n and the WA21 data; 

• the meson-cloud model predicts zero polarization and is in contradiction with all data; 

• the best description of the NOMAD data is achieved in the polarized-intrinsic-strangeness model 
with scenarios Ba and Bb, provided that Dp w 1. The remnant s-quark polarization is higher in 
absolute value for the Aa and Ab scenarios. If one allows a large depolarization during hadroniza- 
tion. Dp ^ 0.4-0.5, then the scenarios Aa and Ab can also provide a fair description of all data. 

It is possible to distinguish between scenarios A and B of the last model by measuring the A polarization 
in the TFR of the neutrino SIDIS. One should expect that P^ ~ i^A in case A and P^ ~ 3Pa in case B. 
The models predict different polarizations for intrinsic s- and s-quark sea: 



57 



NQM: As = As = 0; 



• meson-cloud model []26[]:As « 0, As < 0; 

• intrinsic-strangeness-model A: As « As < 0; 

• intrinsic-strangeness-model B: As w 1/3; As < 0; 

In principle, these predictions can be independently tested by measuring the asymmetries of strange- 
particle production in the current fragmentation region of SIDIS or of (anti)neutrino DIS on a polarized 
target. 

One could try to improve the naive quark model by taking into account the SU (6) symmetry 



breaking as it was done in [168]. In the meson-cloud model one can expect that the contributions from 



higher possible fluctuations with vector meson K~^*K will lead to non-zero A polarization. However, the 



estimates of ref. [ ]26| ] show that the relative probability of this state is small with respect to K^h. (less 



than 10%). Moreover they predict a positively-correlated s-quark spin in a polarized proton. 
9.3 Discussion 

The study of A polarization in the CFR allows us to understand the spin transfer from the quarks to the 
A. Due to the natural helicity and flavour selection of neutrino couplings, we can precisely single out 
specific quark contributions: this information cannot be obtained from the usual lepton-initiated DIS. 
The information on polarized and unpolarized fragmentation functions into A which is available from 
LEP data is scarce and uncertain; also data from NOMAD are far from being decisive in fixing the 
features of quark-spin transfer. The zv-Factory data will induce a big improvement in our understanding 
of fragmentation processes. 

In the target fragmentation region, one can study a new phenomenon, namely, the polarization 
transfer from the lepton to the final hadron. The models discussed here are the first attempts to describe 
this effect. The different models of spin transfer from a polarized quark to a polarized A are able to 
describe the existing LEP data but give different predictions for (anti)neutrino DIS (see, for example. 



[164, 170] and [179] where comparison with NOMAD data are presented). 

Finally, we would like to mention that the best existing data on A polarization, which come from 
the NOMAD experiment, are based on the analysis of about one million DIS events. The statistics at the 
I'-Factory is expected to be a hundred times higher, thus providing ten times better statistical accuracy. 
This will provide a more detailed study of the A polarization dependence on kinematic variables and 
allow better comparisons with different model predictions. 

10 SEARCHES FOR NEW PHYSICS 

We now review briefly some prospects in searching for new physics using the intense beams from a 
neutrino factory. 

10.1 A search for Z' in muon-neutrino-associated charm production 

In many extensions of the SM the presence of an extra neutral boson, Z' , is invoked. A precision study 
of weak NC-exchange processes involving only second-generation fermions is still missing. A search for 



Z' in muon-neutrino-associated charm production has recently been proposed [180]. This process only 
involves Z' couplings with fermions from the second generation. It is interesting because an exotic Z' 
with stronger coupling to the I^, = 1/2 component of weak isospin doublets could still give measurable 
effects at neutrino factories, unlike LHC experiments, which are only sensitive to the Z' coupling to 
charged leptons (/s = —1/2). 

We briefly review the method and the application to neutrino factories using an ideal detector. For 
a detailed discussion of this method, see Ref. [fil 
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10.1.1 The process 

At Q2 ^ Ml the NC effective La grangian ruling the associated charm production induced by and 
including the new physics term is given by 1 180] 



where 



pvcc Gp_ 



ev[c) 



--^i/^7°(l - 75)z^^. c-fa [ey(c) - £^(0)75] c 



€l{u) + eR{u) + 
ev{u) + 



^l{u) - eniu) + 
eA{u) + 



M2 



r]V 



Mf 
Ml 



{rjL + m) 



1 + 



Ml 



Ml, 



VA 



1 + 



'Ml' 



(109) 



(110) 



(111) 



and the parameters x and y give the departure of the couplings from SM predictions. 



10.1.2 Description of the method 

The peculiar topology of the associated charm production in NC interactions is exploited: two 
charmed hadrons in the final state. Consequently, there are no other physical processes that may mimic 
it. Experimentally we are sensitive to the ratio 



(112) 



which can be written as the product 



cr. 



where a^^^ {Z^) is the cross section of the associated charm-production process in interactions in the 
absence of the Z' boson, a^^^ {Z^ + Z') includes the contribution of the new neutral boson, and a^'-' is 
the DIS CC cross section. 



From Eq. (lO) it is clear that the relevant information about the Z' comes from the ratio r, 
which is unity in absence of the Z' . In the following we assume a 50 GeV mono-energetic i/^ beam 
Under this assumption, by using the HERWIG simulation program to compute the the ratio /, we get 
/ = (1.25 ±0.01) X 10~^. 

If we parameterize the ratio r in terms of the x, y and M|, variables defined in Section 10.1.1 , the 
most general expression we obtain is: 

500 _ . / 500 \\ . 2,0 2 



r(x,y,M|,) = l+ (^^j {Aiy + Bix) + ^j^^j {A2y' + B^x' + Cixy). (114) 

Fitting the data from the calculation with the previous function, we obtain the following values of the 
coefficients: Ai = 0.1, A2 = 0.003, Bi = 0.02, B2 = 0.0007 and Ci = -0.0002. The fit is valid in the 
[—30, 30] range for both x and y variables. Figure 32 shows the fitted function r for M^' = 500 GeV/c^. 
The number of observed events, Ns, can be written as 

" The results achievable with a real neutrino spectrum of mean energy {Ei,) are rather well reproduced by using a simple 
mono-energetic beam with energy equal to (-E^). 
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Fig. 32: The ratio r is plotted by assuming Mz' — 500 GeV/c^. 



Ns = N,,-^-r, (115) 

where Ncc is the number of observed events without the Z' effect, es and are the reconstruction 
efficiencies for the events with and without a Z', respectively. 

10.1.3 Measurement accuracy 

Once the charmed particles have been tagged, the Z' effect would show up as an excess/defect of doubly 



charmed events in NC interactions. From Eq. ( |l 14| ) we can argue that for 'large' Z' couplings, i.e. x 



and y > 20, we can get an enhancement of the associated charm production of about a factor 7. On the 
other hand, if we do not observe any excess/defect, we can put a limit on the x and y parameters. As 



an example we report in Fig. |33| the sensitivity plot at 95% C.L. for the x and y variables at Mz' = 500 
GeV/c^. Different systematic errors are assumed from 1% to 50%. 

The allowed region of parameters is obtained from the formula 

1 - 1.96 • <— T < 1 + 1.96- , (116) 

Ncc £B Ncc 



where a is defined as 

^''-stat ^systJ 

and includes the error on the event counting from both a statistical and systematics sources. The factor 
1.96 takes into account the required confidence level. 



In Fig. 33, for each plot, the two lines bound the region of coupling parameters where no significant 
excess/defect of associated charm-production events is found. In other words, an observation of a number 
of charm pair events in agreement with SM predictions excludes the regions outside the band. For each 



plot shown in Fig. we report in Table [llj the Ax and Ay values, respectively the bandwidth at y 
-30 and x = — 30. We assume that 10'' CC interactions are collected at the i/-Factory. Figure 33 shows 



the comparison between the sensitivity achievable with present experiments^ (say CHORUS) and the 

'''Such a meas urem ent could also be exploited by the NuTeV experiment, which recently measured NC charm production in 



-Fe scattering 1 1 8 1 1 
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Scenario 


Present experiments 


i/-Factory 


£ 


Ax 


Ay 


Ax 


Ay 


0.01 


12.0 


5.5 


2.5 


1.5 


0.10 


12.5 


6.0 


3.5 


2.0 


0.25 


14.0 


6.5 


6.5 


3.0 


0.50 


18.0 


8.0 


13.0 


5.5 



Table 11: Band widths at a; = —30 and y — —30 for all the sensitivity plots shown in Fig. e indicates the systematic 
uncertainty. 

one obtainable at a z/-Factory. At a z^-Factory the systematic uncertainties would play an important role 
if they are larger than about 10%, while the statistics is the leading contribution for present experiments. 




Fig. 33: The sensitivity plots for the x and y variables at Mz' = 500 GeV/c^ are shown, e indicates the systematic 
error. 10^ CC interactions are assumed to be collected at a i/-Factory (right plot). The sensitivity achievable with 
present experiments is also considered (left plot). 



10.2 Bounds on 4-fermion operators from a i/-Factory 

Many types of Beyond-the-Standard-Model (BSM) physics appear at energy scales above that of the 
neutrino-scattering process, so their tree-level effects at a i^-Factory can be parametrized by 4-fermion 
operators. These can contribute in the /x^ decay, or in the scattering in the detector. All new physics 
induced by scalars or vectors can be written, via Fierz transformations, as {V ± A){V — A) vertices, 
which we normalize as: 

7/f?f (£7"Pe)(P7a^L^^), Vpfiu^7''PdJ){i^o.PL'^), (118) 

where P = Pl = (1 — 75)/2 or P = Pr = (1 + 75)72, and ^ is a charged lepton. Neutrinos without 
index can be of any flavour. We assume that the strongest neutrino interactions are SM weak, so that the 



neutrino-flavour basis makes sense perturbatively (see ref. [182]). 



Bounds on a variety of extensions of the SM from neutrino scattering were reviewed in [ 183p ; 



a brief update can be found in [142]. New physics in neutrino scattering that cannot be parametrized 
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by 4-fermion operators has recently been discussed in [^]. A catalogue of old constraints on 4-fermion 
operators can be found in Refs. [183] (many models) and [ |184 | (leptoquarks and bileptons) — more up- 



to-date bounds are reviewed in [185] {Z's) and many other recent papers. Constraints on Z's from a 



neutrino factory were studied in the previous section, and bounds on i?-parity-violating couplings from 



neutrino scattering have been discussed in [186] 



A z^-Factory can improve bounds on 4-fermion operators involving neutrinos by many orders of 
magnitude. However, since neutrinos are weakly interacting, the resulting bounds are generically weaker 
than bounds on 4-fermion operators involving charged leptons rather than neutrinos. For instance, new 
physics that can contribute both to lepton-flavour violation in neutrino scattering and to /i — > 67 will be 
more strongly constrained (or more readily detected) in than in neutrino scattering, as can be 

seen from the bounds on lepton-flavour violation discussed in the Report of the Stopped-Muon Working 



Group Q187| ]. It could nonetheless be interesting to look for BSM physics in neutrino scattering at the 
near detector of a neutrino factory. 

First of all, any observed lepton-flavour violation must be due to something other than oscilla- 
tions. Neutrinos are in fact produced only a few hundred metres away from the detector, and have no 
time to oscillate, given the current determinations of the mixing parameters. New physics that can be 
parametrized by 4-fermion operators can then be distinguished from oscillations. 

Secondly, there is BSM physics that can induce uvii 4-fermion operators, without inducing 
operators]^ {£ is a charged lepton.) For quarks, BSM physics (e.g. LQ, Z') that induces uugq also 
induces ££qq, so we do not discuss operators involving quarks. There could nonetheless be interest- 
ing limits on flavour-changing vertices involving charm quarks, from neutrino scattering off nucleons, 
because rare D-meson decays are poorly measured. 

As explained above, we only consider the case of flavour violation induced in the muon decay. 
Since we do not carry out a complete simulation of the event rates, we shall assume that the kinematics 



of the muon decay induced by the first vertex in Eq. (118) is the same as for the SM (V — A){V — A) 
current. We shall then parameterize the muon-decay rate induced by the vertex ( |118[ ) as: 

2^G ) ^ eUeUf,) = [e'^l^'^'fTsMil^ ^ eueiy^) (119) 

and will assume the neutrino beam shapes induced by these BSM decays to be approximated by the SM 
one. 

The new physics signal that we explore is given by a final-state r or by a wrong-sign muon (WSM; 
e.g. a /i"*" produced in the detector by a neutrino beam produced in decay). These signals do not 
suffer from any irreducible physics background. Some r's can be produced in the decays of charmed 
mesons produced by SM CC interactions, but in principle these events will also contain an electron or 
muon from the CC vertex. The r decays from charm produced in NC neutrino interactions, can be vetoed 
by identifying the second charm in the event. Similar considerations apply to events with WSMs. As a 
consequence, a meaningful estimate of the background rates is not possible without a concrete detector 
study. To quantify the discovery potential, we therefore Umit ourselves to presenting a sensitivity reach 
for the detection of N anomalous events, assuming 100% reconstruction efficiency for the BSM signal. 
We assume 10^ ue and 10® uN scattering interactions per year in the detector due to SM CC interactions. 
Then BSM physics could be detected in ue scattering if: 



V 



2V2GF 

and could be detected in uN scattering if: 

V 



2V2Gi 



€ > y/N X 10"^ (120) 



e>\/iVxl0-^ (121) 



'^For example by the exchange of a singlet that couples to {vii 
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Table 12: Flavour changing four fermion vertices involving neutrinos, which a v factory could set a bound on. Neutrinos with 
index i can be of any flavour. The processes listed in parenthesis presently constrain (at 90 % cl) the coefficients of the four 
fermion vertices to be less than e x 2\/2Gf- A neutrino factory could set bounds given in the last column. These bounds 
assume that all 4-fermion operators other than the contrained one are absent. Note that there are much stronger bounds on 
liijik^i vertices than the quoted bounds on liljVkVi- 



Vertex 


Current limits 


i/-Factory limit 


{ll-f°'PLe){Vr-faPLVi) 


e < 0.06 (u oscillations) 


3 X 10-4 




|e| <0.03(G^(g|«)) 


3 X 10-4 


{jl'j'^PLe) {v^-fa Pl v% ) 


e < 0.1 (/u ^ euev'f,) 


3 X 10-4 


{fl-/"PRe){u^-faPLl^i) 


\e\ <0.03(G^(g|«)) 


3 X 10~4 


{ll'y'^PRe){h'aaPLVr) 


|e| <0.03(G^(5|«)) 


3 X 10-4 


ifl^"PLe)iiyaaPLJ^T) 


|e| <0.2(G^) 


3 X 10-4 



Notice that while the reach is lower with ue scattering, the detector backgrounds could be much smaller, 
and any future detailed study will have to consider the possibility offered by this channel. The results we 
present here are obtained assuming the event rates from uN scattering. 

We list in Table [l^ the lepton-flavour-violating 4-fermion operators on which a v factory could 
set better bounds than are available now. In the first column, we list fteuv vertices of different chirality. 
This is a useful way of listing 4-fermion vertices because the new physics that generates them could 
depend on fermion chirality. In the second column are the best available limits yet on the operators, 
and the processes from which the bounds come. In the last column are the limits that a z^-Factory 
could set, assuming a 10 event sensitivity (A^ = 10 in Eqs. 12C and 121). Horizontal lines in the last 
column separate operators that are distinguishable at a i/-Factory because they induce different final- 
state particles in the detector; polarized muons would be required to distinguish between LL and RL 
couplings. Similarly, horizontal lines in the second column separate operators that are distinguishable in 
present experiments (where the flavour of the outgoing neutrinos is often not detected). 



11 CONCLUSIONS 

This work documents our assessment of the physics potential of detectors placed at the front-end of 
a high-current muon storage ring. In most of the cases presented, we tried to evaluate in quantitative 
terms the ultimate accuracies that can be reached, given the available statistics and given the theoretical 
knowledge available today. 

In the case of determinations of the partonic densities of the nucleon, we proved that the z^-Factory 
could significantly improve the already good knowledge we have today. In the unpolarized case, the 
knowledge of the valence distributions would improve by more than one order of magnitude, in the kine- 
matical region x ^ 0.1, which is best accessible with 50 GeV muon beams. The individual components 
of the sea (u, d, s and s), as well as the gluon, would be measured with relative accuracies in the range of 
1-10%, for 0.1 ^ a; ^ 0.6. The high statistics available over a large range of would furthermore al- 
low the accurate determination of higher-twist corrections, strongly reducing the theoretical systematics 
that affect the extraction of as from sum rules and global fits. 

In the case of polarized densities, we stressed the uniqueness of the i/-Factory as a means of 
disentangling quark and antiquark distributions, and their first moments in particular. These can be 
determined at the level of few per cent for up and down, and 10% for the strange, sufficient to distinguish 
between theoretical scenarios, and thus allowing a full understanding of the proton spin structure. A 
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potential ability to pin down the shapes of individual flavour components with accuracies at the level of 
few per cent is limited by the mixing with the polarized gluon. To identify this possible weakness of the 
i/-Factory polarized-target programme, it was crucial to perform our analysis at the NLO; we showed in 
fact that any study based on the LO formalism would have resulted in far too optimistic conclusions. This 
holds true both in the case of determinations based on global fits and on direct extractions using flavour 
tagging in the final state. Our conclusion here is that a full exploitation of the z/-Factory potential for 
polarized measurements of the shapes of individual partonic densities requires an a-priori knowledge of 
the polarized gluon density. It is hoped that the new information expected to arise from the forthcoming 
set of polarized DIS experiments at CERN, DESY and RHIC will suffice. 

The situation is also very bright for measurements of C-even distributions. Here, the first moments 
of singlet, triplet and octet axial charges can be measured with accuracies which are up to one order of 
magnitude better than the current uncertainties. In particular, the improvement in the determination 
of the singlet axial charge would allow a definitive confirmation or refutation of the anomaly scenario 
compared to the 'instanton' or 'skyrmion' scenarios, at least if the theoretical uncertainty originating 
from the small-x extrapolation can be kept under control. The measurement of the octet axial charge 
with a few percent uncertainty will allow a determination of the strange contribution to the proton spin 
better than 10%, and allow stringent tests of models of SU{3) violation when compared to the direct 
determination from hyperon decays. 

The measurement of two fundamental constants of nature, as{Mz) and sin^ 9w, will be possible 
using a variety of techniques. At best the accuracy of these measurements will match or slightly improve 
the accuracy available today, although the measurements at the i/-Factory are subject to different system- 
atics and therefore provide an important consistency check of current data. In the case of as{Mz), the 
dependence of the results on the modehng of higher-twist corrections both in the structure function fits 
and in the GLS sum rule is significantly reduced relative to current measurements, as mentioned above. 
In the case of sin^ 6w, its determination via ue scattering at the i^-Factory has an uncertainty of approxi- 
mately 2 X 10~"^, dominated by the statistics and the luminosity measurement. This error is comparable to 
what already known today from EW measurements in Z° decays. Compared to these, however, this de- 
termination would improve current low-energy extractions, and be subject to totally different systematic 
uncertainties. It would also be sensitive to different classes of new-physics contributions. The extrapo- 
lation to Q = Mz is affected, at the same level of uncertainty, by the theoretical assumptions used in 
the evaluation of the hadronic-loop corrections to ^-Z mixing. The determination via DIS, on the other 
hand, is Umited by the uncertainties on the heavy-flavour parton densities. As shown earlier, these should 
be significantly reduced using the z/-Factory data themselves. 

In several other areas, the data from the i/-Factory will allow quantitative studies to be made 
of phenomena that, so far have only been explored at a mostly qualitative level. This is the case of 
the exclusive production of charmed mesons and baryons (leading to very large samples, suitable for 
precise extractions of branching ratios and decay constants), of the study of spin-transfer phenomena, 
and of the study of nuclear effects in DIS. While nuclear effects could be bypassed at the i^-Factory 
by using hydrogen targets directly, the flavour separation of partonic densities will require using also 
targets containing neutrons. This calls for an accurate understanding of nuclear effects. The ability to 
run with both H and heavier targets wiU in turn provide rich data sets useful for quantitative studies of 
nuclear models. The study of A polarization both in the target and in the fragmentation regions, will 
help clarifying the intriguing problem of spin transfer. We reviewed several of the existing models, 
and indicated how semi-inclusive neutrino DIS wiU allow the identification of the right ones, as well as 
providing input for the measurement of polarized fragmentation functions. 

Finally, we presented some cases of exploration for physics beyond the SM using the z/-Factory 
data. Although the neutrino beam energies considered in our work are well below any reasonable thresh- 
old for new physics, the large statistics makes it possible to search for manifestations of virtual effects. 
The exchange of new gauge bosons decoupled from the first generation of quarks and leptons can be 
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seen via enhancements of the inclusive charm production rate, with a sensitivity well beyond the present 
limits. Rare lepton-flavour-violating decays of muons in the ring could be tagged in the DIS final states 
through the detection of wrong-sign electrons and muons, or of prompt taus. Once again, the sensitivity 
at the i/-Factory goes well beyond existing limits. 

The work presented here has two clear weaknesses, which point to two directions for further work. 
On one side, a realistic evaluation of the experimental feasibility of the proposed measurements should 
be performed. Concrete detector designs should be proposed, and the detector performance should be 
evaluated in the light of the statistical and theoretical accuracy reach set by our study. As part of this work, 
an optimization of the beam parameters (energy, length of the straight section, distance of the detector 
from the ring) should be performed for each individual physics task. On the other side, we made no any 
effort to present a physics case justifying a set of goals for the performances. Considering that these 
results will not be available before at least 10-15 years from now, some judgement on the merit of these 
measurements has to be given. For example, why should we want to know 15 years from now the strange 
density of the proton with an accuracy of 1%? How would our knowledge of fundamental physics, or our 
ability to predict new phenomena, improve if we could reach this goal? From the measurements Usted in 
this document, which ones will be the most important at the time when the z/-Factory will be operating? 
We hope that future studies of this new and fascinating facility will address this important aspect of the 
case for front-end experiments. 
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